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Parafermions are the simplest generalizations of Majorana fermions that realize topological or¬ 
der. We propose a less restrictive notion of topological order in ID open chains, which generalizes 
the seminal work by Fendley [J. Stat. Mech., P11020 (2012)]. The first essential property is that 
the groundstates are mutually indistinguishable by local, symmetric probes, and the second is a 
generalized notion of zero edge modes which cyclically permute the groundstates. These two prop¬ 
erties are shown to be topologically robust, and applicable to a wider family of topologically-ordered 
Hamiltonians than has been previously considered. An an application of these edge modes, we for¬ 
mulate a new notion of twisted boundary conditions on a closed chain, which guarantees that the 
closed-chain groundstate is topological, i.e., it originates from the topological manifold of degen¬ 
erate states on the open chain. Finally, we generalize these ideas to describe symmetry-breaking 
phases with a parafermionie order parameter. These exotic phases are condensates of parafermion 
multiplets, which generalizes Cooper pairing in superconductors. The stability of these condensates 
are investigated on both open and closed chains. 


I. INTRODUCTION 

Parafermions are the simplest generalizations of 
Majorana fermions^ that realize topological order.^ One 
signature of topological order is the manifestation of 
non-Abelian anyons^”^ on a ID open chain,^ which has 
applications in topological quantum computation.^ ^ 
For practical computation, it is desirable that the 
anyons robustly survive any generic local perturbation. 
Following a seminal work by Fendley,^ there is strong 
evidence that topological order is indeed stable over 
a family of ZTv-symmetric Hamiltonians with broken 
parity and time-reversal symmetries. On an open chain, 
these so-called chiral Hamiltonians^^”^^ satisfy an addi¬ 
tional symmetry that renders their entire spectra A'-fold 
degenerate.^ We propose that this symmetry is sufficient 
but not necessary for a topologically ordered phase, since 
topological order describes the low-energy subspace, 
rather than the entire range of excitations. By relaxing 
this symmetry constraint, we are led to consider a 
wider family of topologically-ordered Hamiltonians than 
has been considered in Ref. 2; these include the non- 
chiral Hamiltonian which is dual to the widely-studied 
ferromagnetic clock model.Despite the loss of 
degeneracy in the high-energy subspace, the groundstate 
space remains A^-fold degenerate to superpolynomial 
accuracy in the system size. We are able to demonstrate 
a stronger statement, that the groundstates are mutually 
indistinguishable by all local, symmetry-preserving op¬ 
erators, of which the Hamiltonian is but one case in point. 

The nontrivial groundstate degeneracy on an open 
chain is intimately connected to the presence of edge 
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FIG. 1. We formulate a new notion of twisted boundary con¬ 
ditions on the closed chain. A topological groundstate is espe¬ 
cially sensitive to a change in boundary conditions (schemati¬ 
cally as (a) to (b)). Such a twist results in a different ground- 
state with a distinct Zat charge. 

operators which permute the various groundstates. The 
best-known examples include the dangling Majorana 
fermion on the edge of the Kitaev chain,^^8 
as various other proposals to realize Majorana bound 
states.A related signature is the nontrivial degen¬ 
eracy in the entanglement spectrum, where a real-space 
bipartition introduces a virtual edge.^^’^^ To precisely 
formulate this ‘dangling parafermion’, we introduce the 
notion of a pair of operators (O/,©^), which localize 
on left and right edges of an open chain. Each of the 
pair permutes the groundstates cyclically, but does not 
necessarily commute with the full Hamiltonian, thus 
differing from past proposals.Our generalized edge 
modes are found to generate a non-commutative algebra 
within the groundstate space, which for A" = 2 is the 
well-known Clifford algebra. 

The topological degeneracy that we attribute to 
‘dangling parafermions’ is generically lost on a closed 
chain without edges.One might therefore ask: is the 
non-degenerate, closed-chain groundstate continuously 
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connected to the open-chain topological manifold? In 
this paper, we set out to answer a related question: 
what is the closed-chain Hamiltonian (1-L) that is exactly 
minimized by a state in the open-chain manifold? We 
find that 1-L is obtained by utilizing Q r\j 0\0^ as an 
inter-edge coupling to close the chain. Since Q breaks 
translational-invariance, our proposed Hamiltonian 
differs from conventional methods^^’^^ to close the chain. 
Q may be interpreted as a topological order parame¬ 
ter, in analogy with the local order parameter (O) of 
traditional broken-symmetry phases. Traditionally, by 
applying a local, symmetry-breaking ‘field’ AO + h.c., 
one picks out a broken-symmetry groundstate that 
depends on the phase factor A. Analogously, by applying 
AQ + h.c. to a topological manifold on an open chain, 
one selects a closed-chain groundstate that originates 
from this manifold. By exploring the space of A, we 
are able to select a different state from the topological 
manifold, with a different charge; this variation of 
A is analogous to twisting the boundary conditions in 
integer quantum Hall systems. 


construction of a topological order parameter (Q), which 
lends us a new notion of twisted boundary conditions on a 
closed chain. This notion is formulated in Sec. V; a crit¬ 
ical comparison is made with the conventional method 
of twisting. In Sec. VI, we generalize these concepts to 
exotic phases with parafermionic order parameters. The 
stability of these phases is investigated on both open and 
closed chains. In the last Sec. VH, we discuss the possible 
generalizations of our work. 


II. REVIEW OF PARAFERMIONIC PHASES 
WITH TOPOLOGICAL ORDER 
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FIG. 2. Parafermion chain of L sites. Each site (blue circle) 
is split into two subcells. 


While topologically-ordered parafermions naturally 
generalize the Z 2 Kitaev chain,^ some parafermionic 
phases have no Major ana analog. Indeed, Z^ 
parafermions with V > 4 can also exhibit symmetry¬ 
breaking, sometimes even in conjunction with topological 
order.Symmetry-breaking due to a parafermionic 
order parameter (V) is nonlocal, and thus fundamentally 
different from traditional symmetry-breaking by a local 
order parameter. Expanding on an interpretation in 
Ref. 28, we identify these exotic phases as condensates 
of ‘Cooper multiplets’, which generalize Cooper pairing 
in superconductors. A major advance of Ref. 30 is the 
construction of a class of frustration-free Hamiltonians 
^(iv,n )5 which realize all possible distinct phases. For 
these fine-tuned Hamiltonian, the exact form of y is 
now known (when there is symmetry-breaking), as is the 
exact form of the zero edge modes (when there is topo¬ 
logical order). On the other hand, the stability of these 
phases beyond the frustration-free limit ii?( 7 v,n) has been 
posed as a conjecture.This lack of quantitative under¬ 
standing has led to controversy, e.g., with regard to the 
correct groundstate degeneracy on a closed chain in the 
presence of symmetry-breaking.One goal of this pa¬ 
per is to settle this controversy and prove this conjecture. 

The outline of our paper: we review the simplest mod¬ 
els that realize topological order in Sec. H, and describe 
the strict notion of a zero edge mode on an open chain. 
After arguing that this strict notion is sufficient but not 
necessary for topological order, we then formulate two es¬ 
sential properties of topological phases on an open chain. 
As described in Sec. HI, the first is that the groundstate 
space is mutually indistinguishable to symmetric, local 
probes. The second property is a generalized notion of 
zero edge modes, which we show to have several inter¬ 
esting applications in Sec. IV. One application is the 


Parafermions generalize the Majorana algebra through 

'yf = 71=%^”^ and 7^,7, 7^, (1) 

for / 7 ^ j, and uj = Here, we have denoted 7 ^ 

as an operator acting on subcell j G {1,2,..., 2L} on 
a chain of L sites, in the Hilbert space (C^)®^. Our 
convention is to pair subcells 2j — 1 and 2j into a single 
site jf, as illustrated in Fig. 2. The simplest, fixed-point^^ 
Hamiltonians having topological order are dimerized as 

L-l AT-l 

( 2 ) 

i=i /3=o 

where each ‘even’ parafermion (on an even subcell) is 
coupled only to one ‘odd’ parafermion on the adjacent 
site; these couplings are schematically drawn as black 
lines in Fig. 2 . On an open chain, there is then a dangling 
parafermion on each end which is uncoupled from the 
bulk of the chain. For V = 2, Eq. (2) is the well-known 
Kitaev chain with a dangling Majorana mode on each 
end.^ ^(iv,i) has a Z^v symmetry which is generated by 
the string operator 

Q = (3) 

i=i 

i.e., [i 7 ( 7 v,i),Q] = 0 . Q generalizes the fermion parity 
of Majorana systems through = I. We refer to an 
eigenstate of Q with eigenvalue as belonging to the 
charge sector a G {0,1, 2,..., V — 1} = Zat- Each dan¬ 
gling parafermion O G { 71 , 721 .} is a localized unitary 
operator satisfying 

, O] = 0 and OQ = ujQO. 


(4) 
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These symmetry relations imply that the entire spectrum 
is A^-fold degenerate, where each A^-multiplet comprises 
a state in each charge sector. Fendley defines any oper¬ 
ator satisfying Eq. (4) as a zero edge mode.^ We refer 
to and its equivalence class as purely-topological, 

in the sense that its groundstate space has topological 
order, but is not symmetry-broken by a parafermionic 
order parameter; cf. Sec. VI. As remarked before, 
belongs to a special class of fine-tuned Hamilto¬ 
nians that is frustration-free. To clarify, let us decom¬ 
pose into a sum of two-site operators: = 

"^^=0 ‘ By frustration-free, we mean 

that each of is a mutually-commuting projection, 

and is individually minimized by each groundstate. 
has an intuitive interpretation in the clock representa¬ 
tion of parafermion operators, which we now describe. 
The clock representation is obtained by the generalized 
Jordan-Wigner transformation:^^ 

i-i j 

%j-i = n ^ JJ Tfc, (5) 

k=l k=l 

where (jj and Tj act on site j. These clock operators 
satisfy the algebra 


and — oj^'cFj. ( 6 ) 

By this nonlocal transformation, we find that Eq. (2) is 
dual to the ferromagnetic clock model: 

L-l N-1 

E(vd+i)'’’ (7) 

3 = 1 /3=0 

which has a Z^v symmetry generated by Q = nf=iE - 
To be transparent, let us define an V-dimensional basis 
on each site j, satisfying 

crj\oi) . = uj^\a) = |<^ + 1 mod ( 8 ) 

for a G Zat. Each of projects to a ferromag¬ 

netic alignment of the clock variables on adjacent sites: 
\a)j 0 |<^)j+i- Thus, the groundstate of exhibits 

long-range order with respect to the order parameter 
CTj, i.e., {(TiCr]) is nonzero for large \i — j\. This type 
of symmetry-breaking by a clock order parameter has a 
long history.In this paper, we also describe a differ¬ 
ent type of symmetry-breaking by a parafermionic order 
parameter, which is the subject of Sec. VI. 



f/j 

FIG. 3. Numerical simulation of Z 3 parafermions on a 14-site 
chain. We plot the absolute spread in energy between three 
states {|Tc)}, where is the groundstate of the open-chain 
Hamiltonian (9) in the charge sector o G Z3. 

When deviating from the frustration-free iJ( 7 v,i )5 h is 
not obvious if such localized operators (O) satisfying Eq. 
(4) still exist. Eor illustration, we consider a well-studied 
deformation with Z3 symmetry: 

L-l L 

H = - ^ - / E + h.c., (9) 

i=i i=i 

which reduces to iJ( 3 ,i) (up to a constant) if we set 
f = (j) = 0. The /-term is frustrating, in the sense 
that it cannot be minimized simultaneously with the J- 
term. This frustration tends to destabilize the ordered 
phase, such that localized operators satisfying Eq. (4) no 
longer exist where 0 = 0. However, Eendley has nontriv- 
ially demonstrated that such edge modes persist if 0 is 
deformed away from 0 ;^ this is also understood through 
domain-wall dynamics in the dual clock representation.^^ 
Since both time-reversal and spatial-inversion symme¬ 
tries are broken when 0 7 ^ 0 mod 7 r/ 3 , the resulting 
Hamiltonians are called chiral;it has been suggested 
that breaking these symmetries changes the type of topo¬ 
logical order, or eliminates it altogether.^ This is peculiar 
because the non-chiral Hamiltonian (2) at 0 = 0 exhibits 
the largest spectral gap as a function of 0 , thus one ex¬ 
pects here that the topological properties are most ro¬ 
bust. To substantiate this claim, we numerically evaluate 
the groundstate degeneracy in Eig. 3; this degeneracy is 
shown to be most persistent along the line parametrized 
by (/, 0 = 0). Our point of view is that topological 
order is a property of the groundstate space, and any 
constraint on the excited states is not necessary, Eq. (4) 
being a case in point. This motivates us to formulate two 
essential characteristics of the groundstate space, which 
applies to a wider family of Hamiltonians than has been 
previously considered: (i) a persistent groundstate de¬ 
generacy, despite the loss of degeneracy in the excited 
spectrum (elaborated in Sec. HI), and (ii) the existence 
of generalized zero edge modes which do not satisfy Eq. 
(4), but nonetheless generate a non-commutative algebra 
in the groundstate space, as we will show in Sec. IV. 
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III. LOCAL INDISTINGUISHABILITY IN 
PURELY-TOPOLOGICAL PHASES 

Consider the Z 3 , non-chiral Hamiltonian in Eq. (9) 
with 0 = 0. With / = 0 as well, the entire spectrum 
is three-fold degenerate. With nonzero /, it is known 
that the groundstate splitting is exponentially-small in 
the system size L, while the splittings of the excited 
triplets scales as a power law.^^ This suggests a gen¬ 
eral property of topologically-ordered systems, that a 
nontrivial groundstate degeneracy can persist without 
degeneracy in the excited states. In this section, we 
confirm this general property, and also formulate a 
stronger statement: that the groundstates are mutually 
indistinguishable by any quasilocal, symmetric operator 
(O). By quasilocal, we mean that each term in O has a 
support that is finite-ranged or decaying faster than any 
power, i.e., superpolynomially; see App. A for a precise 
definition of quasilocality. 

For open-chain, quasilocal Hamiltonians that are 
topologically-ordered (without symmetry-breaking), we 
formulate a notion of local indistinguishability in their 
groundstate spaces. We define as the nondegener¬ 
ate groundstate projection in the charge sector labelled 
by a G Ztv, i.e., the groundstate has eigenvalue under 
operation by Q. Let P^n,i) = ^a=o and let O de¬ 
note any quasilocal, symmetry-preserving operator, i.e., 
[0,(5] =0. In the thermodynamic limit, 

P{N,1) O P{N,1) — P{N,1) 1 (10) 

for some complex number c, as we prove in the remainder 
of this Section. Since the open-chain Hamiltonian 
satisfies the same conditions as O, Eq. (10) also implies 
that the groundstate is A-fold degenerate with energy 
c{Hj. 

As a first step, we will prove that Eq. (10) is satisfied by 
the frustration-free P(Ar,i), and then we will extend our re¬ 
sult to frustrated Hamiltonians. Indeed, for finite-ranged 
O, Eq. (10) is satisfied by P(Ar,i) exactly, without finite- 
size corrections. This is most conveniently proven in the 
clock representation (7) of Each of N ground- 

states is a classical product state that is fully-polarized 
in the quantum numbers of aj: \ipcx,i) — \^)j 

a e Zn- The property (10) has a simple explanation in 
the clock representation: local operators have difficulty 
transforming one polarized state 10) into another. To 
substantiate this intuition, we first define an alternative 
groundstate basis which diagonalizes Q. For a G Z^v, let 

!</>«,1) = 7^ (11) 

/3=0 

such that (5|0a,i) = ^^|0a,i); we refer to a as the Z^v 
charge of |0a,i)- Eq. (10) amounts to demonstrating 
that the matrix = (0a,i | | 0 / 3 ,i) is proportional 


to the identity. [O^Q] =0 implies that the off-diagonal 
elements vanish. The difference in diagonal elements, 
“*^^5 reduces to a sum of terms like (0^,i| O 
for // 7 ^ z/. If O has finite range, this quantity vanishes 
because O cannot change the polarization on every site. 
A variant of this proof was first presented in Ref. 38 
to explain the robust groundstate degeneracy of the 
Ising model. If O decays superpolynomially instead, 
has a finite-size correction that we bound 

in App. B. 

While deviating from introduces quantum fluc¬ 

tuations to the classical states {\^cx,i)}^ these fluctuations 
have a length scale that is suppressed by a spectral gap, 
thus preserving the property ( 10 ) to superpolynomial ac¬ 
curacy in the system size. We are assuming a gapped, 
quasilocal interpolation Hg between Hq = and a 

deformed, but topologically equivalent, Hamiltonian Hi; 
the Ztv symmetry is preserved throughout s G [0,1]. By 
gapped, we mean that the spectral gap above the lowest 
N states remains finite throughout the interpolation - 
this implies that the projection Pg to the lowest N states 
is uniquely defined at each s. At this point we do not 
assume that these N states remain degenerate. These 
conditions allow us to define a locality- and symmetry¬ 
preserving unitary transformation Vg , that maps the low- 
energy subspaces as 

Ps = VsPoVl with Po = [V«, Q] = 0. ( 12 ) 

By locality-preserving, we mean that if O is quasilocal, 
so is its transformed version VlOVg. Vg is known as an 
exact quasi-adiabatic continuation;^^”^^ its properties 
are elaborated in App. A, and its explicit form is shown 
in the next Section. Under these conditions, we find 
that if Pq is locally indistinguishable, so is Pg. To prove 
this, let O be quasilocal and symmetry-preserving, then 
PgOPg = VgPoVlOVgPoVl Since V^OVg is also quasilo¬ 
cal and symmetry-preserving, we may apply Eq. (10) 
that we have shown to be valid for the frustration-free 
Pq. It follows that PgOPg = cPg for some constant c, up 
to superpolynomially-small finite-size corrections that 
we bound in App. B. 

Finally, we point out that local indistinguishability is 
a unifying property of many other topologically-ordered 
groundstates, as we elaborate in Sec. VH. 

IV. GENERALIZED ZERO EDGE MODES FOR 
PURELY-TOPOLOGICAL PHASES 

The V-fold degeneracy on an open chain arises solely 
from degrees of freedom on the edges. To crystal¬ 
lize this notion, we would like to construct localized 
edge operators which permute the groundstates cycli¬ 
cally. We already know such a set of operators for 
the frustration-free Hq = P(iv,i)- from (5) and (11), it 
is simple to deduce that ( 0 a,i | 7 i | 0 a+imodjv,i) = 1 and 
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(</'c,i| 72 i| 0 c+imodjv,i) = Given this informa- 

tion, we can construct similar operators for any Hamilto¬ 
nian Hs that is connected to Hq by a gapped, symmetry¬ 
preserving interpolation. This will be accomplished by 
quasi-adiabatic continuation, which is known to preserve 
the matrix elements of operators within the groundstate 
space.In more detail, we decompose the groundstate 
space of Hs as Pg = I? where a denotes 

the Ztv charge of Since [Vs,(5] = 0 for the quasi- 

adiabatic continuation of Eq. (12), the charge of each 
state is preserved. It follows that 

(^a,l|^^|^a+lmodAr,l) = 

{<f>i,l\Or\4>UlmodN,l) = (13) 

with the dressed operators 

Oi=Vs7M and a = V,72.Vt. (14) 

That Oi is localized around site 1 follows from being 
locality-preserving. Since both 71 and are unitary, 
(13) also implies O, 

To gain some intuition about these edge modes, we 
construct Oi for the simplest nontrivial example: a frus¬ 
trated Z 2 Kitaev chain, as modelled by 

L-l L 

Hs=iYl l2jl2j+l + *5 E T'2j-l72j- (15) 

i=i i=i 

Decomposing Hg = Hq + sV^ we point out that Hq is 
the frustration-free Kitaev model. The V term tends to 
destabilize the topological phase, resulting in a mono¬ 
tonic decrease of the spectral gap T^ above the two de¬ 
generate groundstates; it is known that T^ = 2 in the 
frustration-free limit, and T^ = 0, signalling a phase 
transition.^ The quasi-adiabatic continuation can thus be 
carried out for s lying in the real interval [0, 5 ], with 5 < 1 
and > 0. The quasi-adiabatic continuation operator 
has the form of a path-ordered evolution^^ 


we obtain the matrix elements of Vg as: 

i(k;sjV^jj;s} = = (A:;s|as|j;s) (19) 

between states separated by the spectral gap T^, i.e., for 
— Ejl > Fg. This leads to dgPs = i[Ds-,Ps\ and Eq. 
(12). To first order in 5 , the edge mode has the form 

Oi = V,7iVl = 7i + is[Va, 7i] + • • • • (20) 

Employing [i^o, 7 i] = 0 and [V, 72 ^- 1 ] = —2^7^, we are 
thus led to evaluate 




i[Po,7i] = -2i / (itF(2t)e* °^ 72 e 


where we have identified T^ = 2 as the spectral gap in 
the frustration-free Kitaev chain. Now apply that 72 = 
{v_^ -h '^-)/ 2 , where = 72 ± ijs are constants under 
evolution by the dimerized Hq: [Hq^v^] = We 

then find 

/ oo 

dtF{2t) = 73 . 

-CX) 

In the last step, we employed E(±l) = as follows 
from Eq. (18). We thus derive that the edge mode 
spreads as Oi = 7 i + 573 -p ..., and a simple computation 
shows [i^s, 7 i + 573 ] = 0{s^). This result is suggestive 
of a two-fold degeneracy in the entire spectrum, as has 
been derived alternatively in Ref. 2 for higher orders in 
s. However, this property is not generic, as evidenced by 
our next illustration. In App. C, we work out the edge 
mode for the Z3 non-chiral parafermion, as modelled by 
H in Eq. (9) with 0 = 0. We identify the J-term as i^o, 
and /-term as the deformation V. The result is 

= T'l + - w7i71) 

+ 7^*772 73 + 0 (^), ( 22 ) 


Vs = Texp[i f Vg'ds'] (16) 

Jo 

in the Time’ variable 5 , generated by the Hamiltonian 

/ oo 

dt F(rH) (dsH,) (17) 

-OO 

We refer to F{t) as a filter function, whose purpose is 
to cut off the time-evolution of dgHs for large |t|. It 
is desirable that F(t) has the fastest decay for large |t|, 
such that Vg is maximally quasilocal;^^’^^ this is elabo¬ 
rated in App. A. In addition, F{t) is imaginary, so that 
Vg is Hermitian. Let us denote the eigenbasis of Hg by 
Hg\j]s) = By choosing the Eourier transform 

of F{t) as 

/ oo 1 

dtF^^F{t) = -- for \n\ > 1, (18) 


which demonstrably does not commute with Hg. In¬ 
stead, Oi commutes with the projected Hamiltonian 
PgHsPg to superpolynomial accuracy, as follows from 
(i) [ 7 i,Po] = 0 leading to [O^Pg] = 0 , and (ii) the 
groundstate degeneracy shown in Sec. HI. In App. C, 
we have also derived the edge mode Oi for the Z 3 chiral 
parafermion (O< 0 < 7 r/ 3 ). In this parameter range, 
we find to first order in s that Oi may be deformed 
to commute with the full Hamiltonian Hg^ at the cost 
of having a weaker decay away from the edge; this 
deformation involves changing only the filter function in 
the quasi-adiabatic Hamiltonian (17). In this way, we 
recover to first order the ‘zero edge mode’ (Tief^) that 
is alternatively derived in Eq. (32) of Ref. 2 ; we do 
not know if this coincidence persists to higher orders. 
As described in App. D, this deformation of the filter 
function becomes increasingly singular as 0 ^ 0 , i.e., 
the deformed edge mode delocalizes in the non-chiral 
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FIG. 4. Characterizing the XYZ model on a chain of L = 22 sites, with varying parameters Jx and Jy. (a) The spectral 
gap above the two lowest states, for the open-chain Hamiltonian (28). As delineated by gapless lines, the regions B, C and D 
correspond to topological phases, as explained in the main text, (b)-(d) evaluates various types of closed boundary conditions 
by the parity-switch criterion; parameter regions are colored white (resp. black) region where the groundstate parity switches 
(is invariant) upon twisting the boundary conditions. The blue dashed lines are phase transition lines on the open chain, (b) 
Translational-invariant boundary conditions, as described in Eq. (35). (c) First-order topological boundary condition (TBC), 
i.e., Eq. (33) with approximating Eq. (30). This boundary condition is constructed for the topological region C only, (d) 
First-order TBC with defined in App. F2. This boundary condition is constructed for the topological regions B and D 
only. 


limit. In this sense, there is a trade-off between spatial 
localization and commutivity with the full Hamiltonian. 

Let us present another application of these generalized 
edge modes. In the groundstate space, they generate a 
representation of the Zn generator: 

Q = uj(^-^y^OjO^, (23) 

which satisfies = I, and 

may be compared with the representation (3) in the full 
Hilbert space, which operates on each site. The difference 
is that Q is a fractionalized representation, with support 
only near the chain edges.Together, O/, Or and Q 
generate a non-commutative algebra in the groundstate 
space: 

OiOr=ujOrOu OiQ = ujQOu 

and OrQ = o^QOr- (24) 

For N = 2, this is the familiar Clifford algebra, with 
Oi acting as the Pauli matrix in the groundstate 
space, Or as and Q as Clearly, the charge of Oi 
is a topological invariant^^’^^’^^, i.e., it does not change 
under quasi-adiabatic continuation. We say an operator 
O has a charge a if Q^OQ = uj^O. In the present 
context, [Q, Vs] = 0 implies that Oi has unit charge, 
independent of the deformation parameter s. A case 
in point is the Z 3 edge mode ( 22 ), where each term in 
the expansion of Oi has the same charge; the preced¬ 
ing discussion shows that this is true to all powers of // J. 

Finally, the existence of a fractionalized Q implies that 
the groundstate degeneracy is unstable to inter-edge cou¬ 
pling when we close the chain. It is worth refining our 
notion of locality on a closed chain. A case in point is Q 


in the frustration-free limit: 

Q = = o-jcTL H (25) 

i=i 

Since sites 1 and L are adjacent on a closed chain, Q is 
local in the parafermion representation, but nonlocal in 
the clock representation. In short, we call Q parafermion- 
local, but clock-nonlocal. Fractionalization implies the 
existence of a parafermion-local term that distinguishes 
the groundstates. That is, we may legitimately add to the 
open-chain Hamiltonian a term of the form: AQ + /i.c., 
which closes the chain and singles out a unique ground- 
state. Which groundstate is singled out depends on the 
phase factor A, as we elaborate in Sec. V. We remark that 
this refined notion of locality is not necessary for an open- 
chain Hamiltonian, where charge-neutrality imposes that 
all parafermion-local terms are also clock-local; this is 
proven in App. A. Finally, we point out that the local- 
indistinguishability condition ( 10 ) on an open chain ap¬ 
plies only for symmetric, dock-local probes, as may be 
verified in the proof of Sec. Ill; this allows for a clock- 
nonlocal but parafermion-local term, e.g. Q, to close the 
chain and break the groundstate degeneracy. 


V. TOPOLOGICAL ORDER ON A CLOSED 
CHAIN 

The topological degeneracy that we attribute to edge 
modes is generically lost on a closed chain.One 
might therefore ask: what inter-edge coupling guarantees 
that the nondegenerate, closed-chain groundstate (|<L)) 
is topological, i.e., that |^) is continuously connected to 
the topological groundstate manifold of the open chain? 
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Does there exist a coupling that uniquely selects a closed- 
chain groundstate to be exactly one of the open-chain 
groundstates? The goal of Sec. V A is to prove the ex¬ 
istence of this inter-edge coupling for Major ana chains, 
and to describe how the coupling is constructed. In Sec. 
VB, we perform a critical comparison with the conven¬ 
tional boundary condition that is commonly employed. 
We then extend our discussion to Z^v parafermions in 
Sec. VC. 


A. Topological boundary conditions 

The simplest illustration lies in the frustration-free Ma- 
jorana model, whose Hamiltonian on an open chain is 

L-l 

= ( 26 ) 

i=i 

As noted in the previous Section, the groundstate mani¬ 
fold is topologically ordered, and comprises a state | 0 a,i) 
in each charge sector a e Z 2 . We would like to con¬ 
struct a closed-chain Hamiltonian that is exactly 

minimized by the open-chain groundstate This 

is accomplished by identifying the low-energy degree of 
freedom on the open chain, which in this example is the 
complex fermion c = (71 — i 72 i,)/ 2 . It is simple to ver¬ 
ify that | 0 o,i) and differ in their occupation of c. 

Therefore, we are led to = Hq with the 

inter-edge coupling 

(5^7 = i(-l)“7x%. = - 2ctc). (27) 

The phase factor (—I)"" may be interpreted as flux 
insertion, and determines the occupation of the fermion 
(c) in the nondegenerate groundstate. Alternatively 
stated, we obtain two Hamiltonians by imposing periodic 
and antiperiodic boundary conditions for the complex 
fermion Cj = ( 72^-1 — i 72 j)/ 2 , as schematically illustrated 
in Fig. 1. 

Let us generalize our discussion to frustrated models 
that remain topologically ordered. For illustration, we 
may deform the open-chain Hq as: 

L-l 

Hxyz — 23-1123123 + 1123+2 

i=i 

-iJy'l2i-i'y2j+2 + ) • (28) 

We refer to this as the XYZ model, as it can be ex¬ 
pressed as a spin-half model through the Jordan-Wigner 
transformation (5); cf. App. G. As shown in Fig. 4(a) 
for the range | Ja,|, | < 1 (square labeled C), the spec¬ 
tral gap persists above the two lowest-lying states, which 
we denote by in the charge sector a. The low- 

energy degree of freedom on the open chain is directly 
generalized as c = {Oi — iO^j2^ where we have dressed 


the frustration-free edge modes as 71 ^ Oi{Jx^Jy) and 
I 2 L Or{Jx^Jy)^ following Eq. (14). Such a dressing 
exists for any deformation within the square C, where a 
gapped interpolation exists to the frustration-free limit. 
In the presence of interactions, c represents a many-body 
excitation. We then close the chain with the inter-edge 
coupling 

ShU, = i-ir+^Q = (-1)“(1 - 2ctc), (29) 
where Q = —iOiO^ is derived for the XYZ model to be 

Q = 'II 1 I 2 L T 'lJyiHbl2L T 1i12L-a) 

~ Jxilf 2 l 2 . 1 bl 2 L T I1I2L-AI2L-2I2L-1) + • • • , ( 30 ) 

to first order in the deformation parameters; cf. App. F. 
Apparently, Q cannot be derived by translating the bulk 
terms of Eq. (28) to the edge. Since Q represents the 
fermion parity in the open-chain groundstate space (cf. 
Sec. IV), 

= (-i)“+^+'i7:7>, (31) 

for |Ja,|,|J^| < 1. This implies that separately 

minimizes the open-chain Hamiltonian {Hxyz) and 
the inter-edge coupling (SH^yz)^ he., must 

be the nondegenerate groundstate of the closed-chain 
Hamiltonian {Hxyz + ^H^Az)- A corollary is that the 
closed-chain groundstate switches parity on twisting the 
boundary conditions, while the groundstate energy is 
preserved. In short, we call this permutation-by-twisting. 

While we have focused on one model for specificity, 
permutation-by-twisting applies more generally to any 
closed-chain Hamiltonian that may be decomposed as 

=77„,,„ + (-l)“+iQ. (32) 

Here, iLopen is any open-chain Hamiltonian which is 
quasi-adiabatically connected to iLo, and Q = —iOiO^ is 
the fractionalized representation of the fermion parity in 
the groundstate space of i^open- We have thus established 
a bulk-edge correspondence: between the existence of 
topological edge modes on an open chain, to the prop¬ 
erty of permutation-by-twisting on a closed chain. We 
refer to Q as a topological order parameter, in analogy 
with the clock-local order parameter (O) of traditional 
broken-symmetry systems. There, the perturbation 
AO -f h.c. picks out a broken-symmetry groundstate 
which depends on A. The difference is that O breaks the 
symmetry but preserves clock-locality, while Q preserves 
the symmetry and parafermion-locality, but breaks 
clock-locality. 

In practice, the exact form of Q is not easily found. 
However, the first-order truncation (Q^^^) of Q is an¬ 
alytically tractable, e.g., Eq. (30) without the dots is 
an approximation to Q to first order in the deformation 



parameters (Jx^Jy)- We thus describe the closed-chain 
Hamiltonian 

(33) 

as having a first-order topological boundary condition 
(TBC), which we proceed to evaluate. Here and in 
future sections, we use the parity switch as the simplest 
diagnostic to evaluate how ‘topological’ a closed-chain 
Hamiltonian is. Ideally, the nondegenerate groundstate 
parity switches upon twisting the boundary condition, 
for all parameters where the open-chain Hamiltonian is 
topologically-ordered. Indeed, we show in Fig. 4(c) that 
the parity switch occurs (whitened regions) in nearly the 
entire square C. We remark that the boundary condition 
(33) is specifically constructed for the topological region 
C. As we will shortly clarify, regions B and D in Fig. 4 
are also topological, but they require a different set of 
TBC’s. 

To critically evaluate the effect of the first-order dress¬ 
ing, we also consider a boundary condition that is zeroth- 
order in Jx and Jy^ i.e., we apply the inter-edge coupling 
(32) with Q replaced by For this zeroth- 

order TBC, we compute the corresponding parity-switch 
diagram in Fig. 5(a). We find that the parity switch re¬ 
mains robust along — 1 < < 1 , where the open- 

chain groundstates are completely classical; the ground- 
states of the ferromagnetic XXZ model are well-known^^ 
to be fully-polarized. In this classical limit, one may ver¬ 
ify that and Q have identical matrix elements within 
the groundstate space. This coincidence is lost away from 
the XXZ line, where quantum fluctuations play an im¬ 
portant role. Where fluctuations are strongest (in the 
orthogonal direction: —1 < Jx = —Jy < 1), the par¬ 
ity switch does not always occur - the first-order TBC 
significantly outperforms its zeroth-order counterpart. 

1.5 1.5 

1 1 

-To 0 

-1 -1 

-1.5 -1.5 

-1.5 -1 0 1 1.5 -1.5 -1 0 1 1.5 

Jy Jy 

FIG. 5. Evaluation of the zeroth-order topological boundary 
conditions, for the XYZ model on a closed chain of L = 22 
sites, (a) We employ the inter-edge coupling (—I)"" 271721 .. 
(b) Here, we employ (—l)°'Jy 272721 .- 1 . White (resp. black) 
region: the groundstate parity switches (is invariant) upon 
twisting the boundary conditions. 

To conclude our discussion of the XYZ model, we re¬ 
mark that regions B and D (see Fig. 4) also correspond 


to a topological phase, by their quasi-adiabatic connec¬ 
tion to a different dimer Hamiltonian: 

L-l 

^0 ^ V J2j-lJ2j + 2^ 

which has edge modes: 72 and 721 .- 1 . Consequently, the 
TBC in regions B and D differs from Eq. (30), as we 
show in App. F 2 . By the parity-switch criterion, the 
first-order TBC (Fig. 4(d)) is shown to be more effective 
than the zeroth-order TBC’s (Fig. 5(b)). 

B. Comparison with the translational-invariant 
boundary condition 

Having introduced a novel type of boundary condition, 
we would like to make a critical comparison with the 
conventional method to close chains. Assuming that the 
open-chain Hamiltonian (i^open) is translational-invariant 
up to the edge, it is common practice^^’^^ to translate 
the bulk terms of i^open to the edge, and then add phase 
factors for different flux insertions. For the XYZ model, 
this procedure produces the closed-chain Hamiltonian: 

^^XYZ = HxYZ + Jxlil2l2L-ll2L 

- i{-lTJyl2l2L-l + K-^Tlll2L^ (35) 

as we systematically derive in App. E; there, we also 
generalize this procedure for Z^r parafermions. By 
construction, is translational-invariant modulo a 

±1 phase factor. A stronger statement is that 
(resp. is completely translational-invariant in 

the even (resp. odd) charge sector, as we demonstrate 
in App. G. Therefore, we describe as having a 

translational-invariant boundary condition (TIBC). In 
comparison, the TBC generically breaks translational 
invariance, as exemplified in Eq. (30). 

It is commonly believed that if Hopen is topologically- 
ordered, then the closed-chain must also be topolog¬ 
ical, implying that the groundstate of is especially 
sensitive to its boundary conditions.We numerically 
test this expectation for the XYZ model on a finite-size 
chain - see Fig. 4(b). The groundstate parity is invariant 
where quantum fluctuations are dominant, as we have 
also seen qualitatively for the zeroth-order TBC (Fig. 
5(a)); both these boundary conditions are outperformed 
by the first-order TBC (see Fig. 4(c)). A more thorough 
finite-size anaiysis in App. G suggests these conciusions 
are robust in the thermodynamic limit. 

C. Closing the chain for general 'Ln>2 parafermions 

For 'Em >2 parafermions on an open chain, the low- 
energy degree of freedom is no longer a complex fermion 
c. Nevertheless, the N eigenvalues of Q correspond to 
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distinct states in the groundstate manifold, and the gen¬ 
eralized edge modes (O/, Or) induce many-body excita¬ 
tions within this manifold, as per the algebra in Eq. (24). 
To select a topological groundstate on a closed chain, we 
propose the inter-edge coupling: 

1 

= (36) 

/3=0 

which projects to the nondegenerate groundstate with 
charge N — a, for a G Twisting the Hamiltonian 

^ cyclically permutes the charge 

of the nondegenerate groundstate, while preserving the 
groundstate energy. 

To see this, let us label the open-chain ground- 
states as 10 ^ i) in the charge sector a; recall that 
s parametrizes the deformation of the open-chain 
Hamiltonian: Hs = Hq + sV. By construction, 
singles out i) as having the lowest eigenvalue: 

has to be the 

nondegenerate groundstate of H^p + This follows 

from Oi^ Or and Q being unitary, which implies the 
operator norm | | |op < 1 by the triangle inequality. 

Thus, i) separately minimizes H^p and 

For illustration, we return to the Z3, open-chain model 
(9), and now add an inter-edge coupling in the form of 
Eq. (36). In the frustration-free limit (/ = 0), the edge 
modes are 71 and 721 ,, and therefore the topological order 
parameter equals (^* 7 ^^ 721 . • In onr numerical simulation, 
we ignored the dressing of Q for nonzero /, and used its 
frustration-free form throughout; once again, we refer to 
this as zeroth-order TBC. Despite this approximation, 
the zeroth-order TBC correctly identifies the topological 
region in the range of Hamiltonians that we explored; see 
Fig. 6 (a), where the topological region (colored white) in¬ 
dicates a Z3 permutation of the non-degenerate ground- 
state charge, as we twist the boundary conditions. To 
support our claim that this region is topological, we note 
its close overlap with the white region of Fig. 3, where the 
groundstate is three-fold degenerate on an open chain. 

VI. GENERALIZATION TO PHASES WITH 
SYMMETRY-BREAKING 

A new type of phase can occur for Z4 parafermions, 
where the groundstate \0) has a reduced Z2 symmetry. 
This means that \0) no longer has a conserved charge un¬ 
der Q, but it does so with i.e., it is a condensate^^ 
of parafermion pairs. A general classification of Z^v 
parafermions allows for these coherent states, sometimes 
in conjunction with topological order. All distinct phases 
are uniquely labelled by n, which is a divisor of 
The n = 1 phase is purely topological, and has been char¬ 
acterized in the previous Sections; we now turn our atten¬ 
tion to 1 < n < V. It is useful to define m = N/n and the 



0 0.5 1 0 0.5 1 


f/J f/J 

FIG. 6 . Simulation of the Z 3 model (9) on a 14-site chain, 
(a) The chain is closed with an inter-edge coupling = 

+ /i.c., which we call zeroth-order TBC. The 
figure represents the nondegenerate-groundstate charge 
as we twist the closed-ring Hamiltonian through all Z 3 val¬ 
ues of the twist parameter a. White region: Qo is max¬ 
imally permuted, i.e., takes on three values, indi¬ 

cating a topological phase. Gray: two. Black: one. (b) 
For comparison, we have also simulated the translational- 
invariant boundary condition, which involves the coupling 
= - Jcj“-T^^7J72l + h.c.. 

greatest common divisor of m and n: g = gcd(m,n). If 
^ > 1, there is symmetry-breaking with a parafermionic 
order parameter; if m > g, the phase is topological or¬ 
dered; if m > ^ > 1, symmetry-breaking coexists with 
topological order. A major advance of Ref. 30 is the con¬ 
struction of frustration-free models which realizes each of 
these phases on an open chain: 

i=i /3=0 
^ L n—1 

( 37 ) 

jf = l ck:=0 

For example, i^( 4 , 2 ) realizes the above-mentioned conden¬ 
sate of parafermion pairs. i^(jv,n) is a sum of mutually- 
commuting projections, and generalizes the dimer model 
(2). Many groundstate properties are now known in 
the frustration-free limit,including the form of the 
parafermionic order parameters and the topological edge 
modes; these properties are reviewed in Sec. VIA. The 
aim of Sec. VIB is to demonstrate that the phases re¬ 
main stable even if frustrated on an open chain. Sec. 
VIC extends the stability analysis to a closed chain. 

A. Review of general parafermionic phases 

We begin by describing the groundstate properties of 
i^( 7 v,n )5 of which many are known from Ref. 30, then 
extend them to more general Hamiltonians by quasi- 
adiabatic continuation. It is instructive to interpret 
H'( 7 v,n) in the clock representation, where all nontrivial 
phases arise from traditional symmetry-breaking with a 
clock-local order parameter. For illustration, we con¬ 
sider the clock-analog of the parafermion-pair conden- 
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sate. One possible groundstate of i^( 4 , 2 ) is fully-polarized 
as (8)^=i(|0)j + |2)j)- This state breaks the Z 4 symme¬ 
try under \a)j \a 1 mod 4)^, but retains a reduced 
symmetry under \a)j \a -\-2 mod 4)^. More generally, 


H 


{N,n) 


1 

m 




i=l /3=0 


^ L n—1 

08 ) 

j=l q :=0 


Ztv/^, i.e., each (m/^)-multiplet has a residual symmetry 
generated by If m > ^, there exists a remnant 

{m/g)-fo\d degeneracy which originates from topological 
edge modes. These localized operators and 7 ^^) 
commute with and permute the groundstates as 

1^/3, 5 ) ^ 1^/3, 5-1 modm/g)- Ih App. H, wc cxcmplify this 
discussion with the model i^( 8 , 2)5 which has coexisting 
topological and symmetry-breaking orders. 


has m classical groundstates |' 0 a,n )5 which are fully po¬ 
larized in the quantum numbers of 



where a G and CFj\(y). = uj^\a )The Z^r generator 
(3) acts in this basis as QI'0c.,n) = l^c+i mod m,n)5 as follows 
from (8). In the charge eigenbasis, 

^ m—1 

= (40) 

^ /3=0 

with a G Z^ and Q|0c.,n) = Since |0c.,n) 

are fully-polarized in the clock representation, clock- 
local operators have difficulty transforming one |0) into 
another; this transformation would involve turning the 
clocks on all sites. This observation leads to the local- 
indistinguishability of the groundstate space P(^ = 

E:=7I0»,.)(</>.,J, i-e., satisfies the same condi¬ 

tion as P(jv,i) in Eq. (10), when probed by any symmet¬ 
ric, clock-local operator (O). To show this, we first sup¬ 
pose the probe O has a finite range. We would like to 
show that the matrix = (0«,n| | 0 / 3 ,n) is propor¬ 

tional to the identity. The off-diagonal elements van¬ 
ish because [O, Q] = 0. The difference in diagonal ele¬ 
ments, equals a sum of terms proportional 

to ('0^,n| O \'ipu,n)^ with fi ^ V. This quantity vanishes due 
to the above-mentioned observation. If O is not finite- 

ranged but decays superpolynomially, van¬ 

ishes up to finite-size corrections that we bound in App. 
B. An instructive basis for P^N,n) is 

^ (41) 

V ^ a=0 

where amjg S E Z^, G Z^ may be interpreted as 
a broken-symmetry index, and 6 G ’^m/g ^s a topolog¬ 
ical index; recall g = gcd(m,n). On application of a 
symmetry-reducing on-site ‘field’: + h.c., the m- 

multiplet splits into g number of (m/^)-multiplets, as 
determined by 

727-i|%. 5) = and 

7277/3,5) (42) 

For any subcell j, are parafermionic order pa¬ 

rameters which reduce the Z^v symmetry of P(jv,n) to 


The exact form of \0) allows us to expand upon their 
interpretation^^ as coherent states. Eq. (42) implies that 
1^) is a condensate of the operator (pj = 7 /^^. This con¬ 
densate is not of Bose-Einstein type, since 

[ 7 .”'^^®, 77 ^®] =0 for e Z, V i,j, (43) 

implying that = 0 ^ Sij. Never- 

theless, off-diagonal long-range order manifests as 
1^/3, <5) = 1, independent of \i - j\. Note 
that [cpi^cpj] = [cpl^cp]] = 0, as would a bosonic opera¬ 
tor. The present situation is reminiscent of long-range 
order in the BCS wavefunction, where the order param¬ 
eter is almost bosonic. There are differences: (a) the 
groundstate expectation 1 ^/ 3 , <5) takes on g discrete 

values (cf. Eq. (42)), in contradistinction with conven¬ 
tional BCS wavefunctions that have a U{1) degree of 
freedom, (b) The order parameter (p manifests an attrac¬ 
tion between N/g parafermions (^ > 1), thus generalizing 
Cooper pairs to ‘Cooper multiplets’ if N/g > 2. 


B. Robustness of general parafermionic phases on 
an open chain 

Let us address the stability of these phases as we 
symmetrically deform iL(iv,n) to a new Hamiltonian Hg. 
As long as the deformation preserves the gap above the 
lowest m states, there exists a quasi-adiabatic contin¬ 
uation (Vs) which maps their respective groundstate 
spaces as P^N,n) to P^. Thus if P(iv,n) is indistinguishable 
to clock-local probes (as shown in Sec. VIA), then this 
property robustly carries forward to P^, following a 
simple generalization of Sec. III. Moreover, in close 
analogy with Sec. IV, we find for ^ > 1 : generalized 
order parameters Vj = which are dressed 

versions of (pj = Dressing Vj may be interpreted 

as spreading the Cooper-multiplet wavefunction, in 
contrast with the tightly-bound pj. By construction, Vj 
permutes the groundstates of Pg, in the same manner 
that Pj would for While [ 7 ?^ , = 0, Vj 

commutes only with the groundstate-projected Hg. 

If the deformed Hamiltonian Hg is topologically- 
ordered, there exist generalized edge modes {Opn^Or,n} 
which permute the groundstates, and are related by 
quasi-adiabatic continuation to { 7 ^: 7?!.}5 here, the sub¬ 
script I (resp. r) indicates that the operator is localized 
on the left (resp. right) edge. These edge modes generate 
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a fractionalized representation of the generator: 

Qg (44) 

with integer x uniquely satisfying xn = g mod m. This 
operator acts like in the groundstate space of 
i.e., in the basis (41), 

Qg\O0,s)=UJ<>^^\O0,s). (45) 

Decomposing Qg into two operators with support on 

opposite ends of the chain, the charge {—x) of the 
left-localized operator ((Oz,n)~^) is a topological in¬ 
variant. Within each (m/^)-multiplet labelled by (3, 
{Oi^n^ Qg} generate a non-commutative algebra: 

Ol^nOr,n Or,nOl^n, Oi^nQg = QgOl^n, 

and Or,nQg=UJ^^QgOr,n^ (46) 

By our assumption of topological order (m > g), the 

phase factors and are never trivially unity. For 
a detailed discussion of this algebra for N = S and n = 2 , 
we direct the interested reader to App. H. 


C. Robustness of general parafermionic phases on 
a ring 

On a ring, the degeneracies due to topological edge 
modes are generically lost; for ^ > 1 , there remains a g- 
fold degeneracy on the ring due to broken symmetry. The 
instability of the topological degeneracy may be under¬ 
stood in this way: the edge modes themselves imply the 
existence of a Z^v singlet Qg which splits the topological 
degeneracy. Suppose we couple the edges as 

= -uj^^^Qg + h.c., with a G Z^/g. (47) 

This coupling splits the open-chain m-multiplet into 
{m/g) number of ^-multiplets, as indexed by ^ G Z^/^ in 
Eq. (45). Indeed, the m-multiplet is a condensate of both 
topological edge modes and symmetry-breaking order 
parameters, while each ^-multiplet is only a condensate 
of the symmetry-breaking order parameters. By twisting 
the inter-edge coupling as a ^ a + 1 , we permute 

the groundstate multiplet as^ = m/g — a^m/g — a — l. 

It is interesting to determine if this ^-fold degeneracy 
persists as we frustrate the closed-chain Hamilto¬ 
nian. Ref. 28 claims that the degeneracy is broken, 
while Ref. 30 suggests that the degeneracy is robust, 
though without proof. The intuition for a general 
proof can be obtained from the simplest symmetry- 
broken model ( 77(4 2 )) without topological order. Here, 
m = n = g = 2. The Z 4 symmetry is reduced to Z 2 due 
to the parafermion order parameter 7 ^^, which manifests 
in a two-fold-degenerate groundstate (|' 0 c.ez 2 , 2 ) from Eq. 
(39)) on both open and closed chains; we denote this 
groundstate projection by P( 4 , 2 )- We would like to know 


if this degeneracy is stable under deformations of 77 ( 4 , 2 ), 
as would be implied if we prove the indistinguishability 
condition: POP = c{0)P; here, P is a projection 
that is quasi-adiabatically connected to P( 4 , 2 ), and O 
refers either to a perturbation, or to the Hamiltonian 
itself (which shall also be termed a ‘perturbation’ 
in the following discussion). A legal perturbation is 
parafermion-local and charge-neutral; all such perturba¬ 
tions on an open chain are also clock-local, as proven 
in App. A. We have shown in the preceding Section 
that the indistinguishability condition is satisfied with 
clock-local perturbations, thus the degeneracy is stable 
on an open chain. However, a closed chain allows for 
inter-edge perturbations which are parafermion-local but 
clock-nonlocal, a case in point being 7 ^ 721 . = 

What remains is to demonstrate the indistinguishability 
condition under clock-nonlocal perturbations, which we 
shall first carry out in the frustration-free limit. A conve¬ 
nient intermediate step is to express in the charge 

eigenbasis: |0«,2) = (|'0o,2) + jj-l)"|'0i,2)_)/y2. Then 
the off-diagonal elements of = ( 0 ^, 2 1 O | 0 / 3 , 2 ) vanish 
because O is charge-neutral while different |0) carry 
different charge. Our frustration-free proof is complete 
once we demonstrate the vanishing of the difference in 
diagonal elements: = (' 0 o, 2 | O |' 0 i, 2 ) + h.c.. 

Returning to the example of 7 ^ 721 ., while Q transforms 
l^/^i 2 ) into |' 0 o, 2 ), 2 ) is orthogonal to the ground- 

state space. More generally, any odd power of aj (for 
any j) takes us out of the groundstate space. This 
follows from = (-l)"^a/ 3 ^a,omod 2 , wMch 

in turn is deducible from Eq. (39). This motivates us 
to try an inter-edge coupling with an even power of 
O = 7 ^ 72 ^ ^ However, any even power 

of Q acts trivially as ('0a,n| |' 0 / 3 ,n) = thus O 

cannot distinguish between |' 0 i, 2 ) and |' 0 o, 2 )- Einding a 
parafermion-local operator that would do so turns out to 
be impossible, as we now show. Any inter-edge coupling 
may be decomposed into the form: , where 

Bi (resp. Br) is a clock-local operator with definite 
charge, and with support near site 1 (resp. L). The 
Jordan-Wigner transformation (5) implies that there 
are as many powers of the string Q as the charge of 
Br. To transform |' 0 i, 2 ) into |' 0 o, 2 ), one needs an odd 
power of the string operator Q; however, an odd power 
of Q always accompanies an odd power of aj (in Br), 
which takes us out of the groundstate space. This 
intuitive observation in the frustration-free limit can be 
supplemented with quasi-adiabatic techniques to prove 
the indistinguishability condition in the presence of 
frustration, as we show in App. I. 


We now extend our discussion to more general {N, n)- 
phases with a parafermionic order parameter, possibly 
in conjunction with topologically order. Let us de¬ 
compose the groundstate space P[H^p] of an open-chain 
Hamiltonian into eigenstates of the Z^v generator: each 
of ^ couscrvcd charge na. Eor ^ > 1, 
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the parafermionic order parameter (or its dressed 
version yj) permutes the groundstates as ^ 

as evidenced from = uj^^\ 

Thus P[i^op] divides into m/g number of parafermion- 
condensed subspaces, which we label by ^ G 1>rn/g arid 
denote as 

^-1 

^<5 = ^ \(t^arn/g+5,n){^arri/g+5,n\' ( 4 ^) 

Ck:=0 

We propose that each of Ps is indistinguishable under 
any symmetric, parafermion-locsl probe 0 , i.e., in the 
thermodynamic limit, 

PsOPs = c{d)Ps (49) 

for a complex number c. The proof of Eq. (49) generalizes 
our intuition from the P( 4 , 2 ) model, and may be found in 
App. L We note that Eq. (49) is a stronger statement 
than the local-indistinguishability condition ( 10 ), which 
applies only to symmetric, clock-locsl probes of P[H^p]. 
By ‘stronger’ indistinguishability, we mean that a certain 
groundstate space {Ps) cannot be distinguished by both 
clock-local and -nonlocal operators (which are symmet¬ 
ric and parafermion-local). Suppose we now close the 
chain symmetrically; any legal term in the closed-chain 
parafermion Hamiltonian satisfies the same conditions on 
0. Since each of Ps satisfies the strong indistinguishabil¬ 
ity condition, but their direct sum P[H^p] = ^ 

does not, the open-chain groundstate splits into multi- 
plets. A multiplet indexed by 6 is robustly degenerate, 
so long as the gap persists within each of the g charge 
sectors of Ps. 


VII. DISCUSSION 

Topological order is known to characterize the low- 
energy subspace for a variety of condensed-matter 
systems.In this paper, we have precisely char¬ 
acterized the groundstate of ID parafermionic chains, 
and identified the properties which are robustly associ¬ 
ated with topological order and/or symmetry-breaking. 
Our work generalizes a previous work that imposed 
constraints on the high-energy states.^ These constraints 
cannot be realized in nature, since parafermions only 


emerge as low-energy quasiparticles. 

A unifying property of topologically-ordered ground- 
states is their indistinguishability by local probes. Erom 
this fundamental property, one may derive:^^ (i) a 
nontrivial groundstate degeneracy depending on the 
topology of the manifold, and (ii) well-known signatures 
in the entanglement entropy.Besides parafermions, 
a close variant^^’^^ of Eq. (10) applies to other models 
with topological order, including the toric code^ and 
Levin-Wen string-net models;a stronger version of Eq. 
( 10 ) is known to stabilize the spectral gap of frustration- 
free Hamiltonians under generic local perturbations, 
and produces an area-law for the entanglement entropy. 

A few works^^’^^ have demonstrated that Majorana 
edge modes may exist with weaker localization properties 
in particle-number-conserving superconductors, where 
the U{1) symmetry (associated with the electron charge) 
is not spontaneously broken to Z 2 . Eor one particular 
soluble model, it has been shown that the groundstate of 
the topological phase switches parity when the boundary 
condition is changed;it is interesting to determine 
if this property is more generally true away from the 
soluble limit, perhaps with similar techniques that are 
presented in this paper. 

In the final stages of this work, the phase diagram for 
our Z 3 model has been alternatively derived from an 
entanglement perspective in Ref. 54. 
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Appendix A: Local operators and quasi-adiabatic continuation 

To define a local operator (O), we first decompose it as 

= Y E (Ai) 

r>lAeS(r) 
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where Vr^A acts nontrivially on the chain A of length r, and S{r) is the set of all chains with length r. For a chain 
of length L, O acts in the Hilbert space: Z = The operator norm ||>C||op of an operator jC is defined by the 

supremum of {||>C'i;|| : v G Z with \ \v\ \ = 1}. We call O quasilocal if it can be characterized by strength K and decay 
/(r), such that the operator norm ||W,a||op < Kf{r) for all r > 1 and for any chain (A) of length r; / : i-G [0,1] 

and decays superpolynomially in r, i.e., faster than any power. Included are finite-ranged and exponentially-decaying 
interactions. For example, the strength of the finite-ranged Hamiltonian (2) is J. In this example and the rest of the 
paper, we consider only quasilocal Hamiltonians, which are characterized by a finite Lieb-Robinson velocity,^^ i.e., 
information effectively propagates at finite speed, in analogy with the speed of light. This implies that any finite-time 
evolution by a quasilocal Hamiltonian preserves the locality of an operator.We now elaborate on an important 
example of such an evolution: a quasi-adiabatic continuation. For illustration, let us symmetrically deform of 

Eq. (2), while preserving the spectral gap above the lowest N states (which we do not assume to be degenerate). We 
assume there exists a family of symmetry-preserving Hamiltonians Hg, which are differentiable in 5 for 5 G [0,1], and 
Hq = The spectral gap allows us to uniquely define Pg as the projection into the lowest N states of Hg. We 

define an exact quasi-adiabatic continuation Vg that maps Pg = is the unitary transformation (16) with 

the Hermitian Hamiltonian (17). The infinite-time evolution in Eq. (17) is filtered by a function F{t) that decays 
faster than any power, for large |t|; an explicit construction of F{t) is provided in the appendices of Ref. 40 and 39. 
We assume for simplicity that the interactions in Hg decay exponentially, or faster. It follows from Lemma 2 of Ref. 
40 that Vg is a quasilocal Hamiltonian. By interpreting s G [0,1] as a time variable. Lemma I of Ref. 40 informs 
us that the finite-time evolution Vs generated by Vg preserves locality, i.e., V^OVg is a quasilocal if O is quasilocal. 
Since [Hg,Q] = 0, it follows that the quasi-adiabatic continuation is symmetry-preserving, i.e., [Vg,Q] = [Vs,Q] = 0. 


It is useful to distinguish between two notions of locality on a closed chain. We say that an operator is clock-local 
(parafermion-local) if it has strength J and superpolynomial decay f(r) in the clock (parafermion) representation. 
We describe parafermions (clocks) only with a parafermion-local (clock-local) Hamiltonian. We also insist that 
all Hamiltonian terms are charge neutral, by which we mean they commute with the Z^v generator. On an open 
chain, we now demonstrate that all parafermion-local terms are also clock-local, due to the charge-neutral constraint. 
Decomposing a generic parafermion-local term as in Eq. (AI), we consider each Vr^A separately: 


(X 7l 


7 

li+i 


+ 1 

• /i+r—1 ’ 


(A2) 


where rij G Z^v, and by assumption of an open chain, i >1 and i -\- r — 1 < L. Let us express this term in the clock 
representation through Eq. (5). Loosely speaking, this transformation introduces strings of Tj for j < i, assuming 
i > 1. However, if Vr^A is charge neutral, the strings all cancel out. More precisely, we operate on each site labelled 
hy j <i as with n = However, 


Q^Vr,AQ = 


(A3) 


and charge neutrality imposes that [Q, W,a] = 0, or equivalently, n = 0 mod N. 


The discussion thus far applies to an open chain. On a closed chain, charge-neutral inter-edge couplings can be 
either (i) parafermion-local but clock-nonlocal (e.g. Eq. (25)), or (ii) clock-local but parafermion-nonlocal (e.g. crja^). 
Refining our notion of locality thus distinguishes between parafermions on a ring and clocks on a ring. 


Appendix B: Indistinguishability by general quasilocal probes 


In Sec. HI, we claimed that (a) P(n,i) satisfies the local-indistinguishability condition (10), for any Ziv-symmetry- 
preserving, clock-local operator O. More precisely, O has strength K and superpolynomial decay /(r) in the clock 
representation, as we describe in App. A. Later in Sec. VIA, we also claimed that (b) P^N,n) is locally indistinguishable. 
The proofs of both statements can be combined. In what follows, n is a divisor of V, and m = N/n. It suffices to 
show that the matrix = (0a,n|C^|0/3,n) is proportional to the identity, where |0a,n) ^ P[H(N,n)] is defined in Eq. 

(40). The off-diagonal elements vanish because [O, Q] = 0 by assumption. The difference in diagonal elements is 




00 


an(iJ,- 


■y) _ ^/3n(^-zy) 


m 




(Bl) 


We claim that 


^2(m-l)X/(L) 


(B2) 
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for a chain of length L. 

Proof: For O to have nonvanishing matrix elements between ('0^,n|C^|'0zy,n) foi* /i 7 ^ z/, we require that O flips the 
clock variable on all L sites. Thus if we decompose O as in App. A, we find ('0^,n|C^|'0zy,n) = where 

Vl,a acts nontrivially on the entire chain. By the Cauchy-Schwartz inequality, 


^L,A 


< 




VL,A\f^v,n) 


VL,A\f^v,n) 


(B3) 


Now applying the definitions of the operator norm and quasilocality (cf. App. A), 


VL,A\f^v,n) 


< 


Vla 


<Kf{L). 


It follows from this and Eq. (Bl) that 

{0a,n I 1 0Ck:,n) {0/3,n | | 0/3,n ) 


- jyi ^ ('0/x,n|O|'0zy,n) 




<2{m-l)Kf{L). 


(B4) 


(B5) 


Appendix C: Generalized edge mode for the Z3 parafermion model 

Let us evaluate the zero edge mode Oi for the Z 3 model (9), which we rewrite as Hs = Hq sV, with 5 = //J, 

L-l ^ L 

=-A y] ^ 72 ^ 172 ^ +/i.c., A = and V =+ h.c.. (Cl) 

1=1 i=i 

We need certain relations introduced in Sec. IV, including the definitions of the quasi-adiabatic continuation operator 
Vs in Eq. (16), which is generated by the Hamiltonian Vg in Eq. (17). The form of Oi in Eq. (20) motivates us to 
evaluate the operator ^[Po^Ti]- Employing [i7o,7i] = 0 and [V, 7 i] = (cj* — uj){j 2 — ^ 7172 ), we are led to 

i[-Do,j,] = {oj* - oj){X - ojjIY) (C2) 


with ^ = j dtF{3t) 7 ^ and ^ = J dtF{3t) 7 . 




(C3) 


Here, we have applied that the spectral gap T^ (above the lowest three eigenstates of Hq) equals 3. The problem is 
reduced by noting Y = — Vi, as follows from F being imaginary. Eurther progress is made by realizing that 72,73 
and form a closed linear algebra under commutation by Hq. That is, the coefficients defined by 


_ Hq, ^o72 + -B 073 + Co ujAiI ] = MI 2 + Bi 73 + Cl ^771 


are linearly related as 


'AC 

B\ I = 3A1 I Bq 


with Hermitian matrix A4 = 


.Cl, 


.Co. 


iVS 


0 -A A* 
A* 0 -A 
-A A* 0 


Diagonalizing A4, we find eigenoperators: 

2 

Vj = 7(^2 + with corresponding eigenvalues Cj{4>) =-p sin + ^j), 

^ V 3 

such that [Hq^vj] = Expanding 72 in this eigenbasis as 72 = ('Cq Y vi Y V 2 )/V^. 


i=o 


373 


(C4) 


(C5) 


(C 6 ) 


(C7) 


t=i 
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To proceed, we thus need the Fourier components F for the frequencies How are they determined? The two 

essential properties of a quasi-adiabatic continuation is that (a) it evolves the groundstate space Pg as a function 
of the deformation parameter 5, and (b) that it is locality-preserving. The first property (a) determines the Fourier 
components F(Q) = —1/fl for |fl| > 1 (as shown in Sec. IV), while the second (b) determines the components for 
\Q\ < 1, up to a degree of arbitrariness. As will shortly be clarified, this arbitrariness lies in the operators which have 
zero matrix elements within the groundstate space. Now we address separately the chiral and nonchiral cases: the 
nonchiral limit 0 = 0 is discussed in App. C 1, and O<0<7r/3is elaborated in App. C 2. 


1. Nonchiral Z 3 parafermions: 0 = 0 

As defined in Eq. (C6), the eigenvalues of M. are 0o = 0, Ci = — 1 and C2 = 1- For the first eigenvalue, the oddness 
of F implies that F{(q) = 0; for j G {1,2}, F{(j) = —l/Cj since \(j\ > 1. Eq. (C7) leads to 

Inserting the expressions of X and Y = —X^ into Eq. (C2), we thus derive from Eq. (20) the form of O/, to first order 
in the deformation parameter s. This result is found in Eq. (22). 


2 . Chiral Z 3 parafermions: 0 < 0 < 7r/3 


In this range of 0, both |0o| and |0i| < 1, while |02| > 1- To determine the edge mode, we thus need the Eourier 
component F{ft) for 0 < |fl| < 1, which enters the expression (C7). In the conventional quasi-adiabatic continuation, 
F in this range is optimized so that Vg maximally preserves locality, which amounts to finding F{t) with the fastest 
possible decay. This has been proven to be of sub exponential type,^^’^^ that is, F{t) decays faster than exp— 
for any a < 1; an example of F is described in the next paragraph. We refer to subexponentially-decaying edge 
modes {Oi) as maximally localized. On the other hand, one may choose an edge mode Oi with worse localization 
properties, but with the desirable property that it commutes with the Hamiltonian. There is thus a trade-off between 
localizability and commutativity. Indeed, choosing F{Q) = —1/Q for \Q\ = \Q\ < 1 worsens the decay of F{t) (and 
thus of the resultant edge mode), but leads to an edge mode which commutes with the Hamiltonian, as we elaborate 
in App. D. One may verify that the difference A = Oi — Oi has zero matrix elements in the groundstate space. To 
illustrate this for the Z3 chiral example, we have shown that the first-order correction to the edge mode has a term 
proportional to T(Ci)'^j; this follows from Eq. (C2) and (C7). Since |0o|, |Ci| < 1? the coefficients of vq^vi (as 

defined in Eq. (C6)) indeterminately depend on our choice of the filter function in the range 0 < |fl| < I; once again, 
we are choosing between spatial localization or commutivity with the full Hamiltonian. This lack of determinacy has 
no physical consequence in the groundstate space, as Vj\2pa,i) = 0 for j G {0,1} and a G Z3; here, recall 0^ are the 
groundstates of i7o, as defined in Sec. HI. 


To aid future numerical studies, we provide further details of the optimized filter function. Here, we do not claim 
originality, but hope to present an accessible summary of certain portions in Ref. 39 and 56. Expressing 


m) 


i-m) 

n 


(C9) 


we define a real, even function g to satisfy: (i) ^(0) = I, and (ii) g has compact support within the interval [—1,1]. 
The uncertainty principle then states that g cannot be ‘very small’ - optimally,^^ g{t) < C exp[—|t|5(|t|)] for some 
constant C, and 5 a positive, monotone-decreasing function satisfying that 


e(y) J 

- ay converges. 

y 


It is further shown by Ingham in Ref. 56 how to construct g such that 

e{t) = 


(CIO) 


(log(2 + i))2 ’ 


(Cll) 
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though it has been noted that faster asymptotic decays with different 5 are possible in principle.For Ingham’s 
choice, 


with 


9{t )=^ n 

n=l 


sin (pnt) 

Pn^ 


(C12) 


Pn — 


Be^np) 

no ’ 
Bs{n) 
n ’ 


n <no 
n > no. 


(C13) 


^4, B and no are chosen such that g satisfies the above properties (i) and (ii); the final result is shown in blue in Fig. 

7. 



FIG. 7. Blue: optimized filter function (F(fl)) obtained by roughly estimating the parameters in Eq. (C13) as A = 
1/9.2843, B = 80, no = 100. Red: — 1/Q is shown for comparison. 


Appendix D: Zero edge modes that commute with the Hamiltonian 

Let us define the first-order truncation of Oi as = 7i + where s is the deformation parameter. This 

operator depends on our choice of the filter function F, which enters the expression for Vg in Eq. (17). For a specific 
choice of F (with corresponding edge mode we find that and Hg commute to first order in s. The following 

discussion clarifies this choice of F. Let Vj be operators which satisfy two properties: (i) they belong in the charge 
sector Of, i.e., Q^VjQ = and (ii) they are constants under time evolution by iLo, i-e., [Ho^Vj] = 2(jVj with real 

eigenvalues 2(j. These operators vj span the space (T/^) of operators in the same charge sector, i.e., for each element 
X G Tb, Q^xQ = uj^x. Labelling each element by a subscript (e.g., x^), commutation by iLo may be interpreted 
as an eigenvalue problem: [Ho^xi] = BijXj. The matrix B is known to be Hermitian in the frustration-free limit,^ 
where Ho = Fat,!, as described in Eq. (2). Applying the completeness property, we expand [V^,71] = with 

c-numbers r^. This leads to 

*[^^0,71] = F ^ ^F( 2 C„/ro). (Dl) 

n n ^ 

If the Fourier transform F(Q) is chosen as —l/Vt for all relevant frequencies Q. G {2C^/ro} in the above sum, we 
denote the resultant edge mode by 0/. This choice of F, which we call Fendley’s choice, is only possible when none 
of the relevant frequencies are vanishing, i.e., Cn 7^ 0 for any nonzero The proof of commutivity to first order is 
elementary: 

[H,,dP] = s[y,7i] + s[Fo,i[Po,7i]] + 0(s2) = s[l/,7i] + s V ^d:^l^F(2C„/ro) + 0{s^) 

n ^0 

= .[y, 71] +. E ( - ^) + 


(D2) 
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In cases where all the relevant frequencies are large, by which we mean they satisfy: |2C^/ro| > 1, Fendley’s choice 
is not so much a choice but a necessity. Indeed, Fendley’s choice then coincides with Eq. (18) in Sec. IV; we have 
shown therein that Eq. (18) leads to the correct evolution of the groundstate space: dsPs = Ps]- The set of cases 
where all the relevant frequencies are large include all quadratic Majorana models, as we show in App. D 1. On the 
other hand, if some of the relevant frequencies satisfy: \2(n/ro\ < 1, Eendley’s choice produces an edge mode that 
is deformed from the maximally-localized edge mode. In App. D 2, we carry out this deformation for the Z3 chiral 
parafermion. 


1. Zero edge mode for quadratic Majorana models 

Suppose we deform the topological dimer model described in Eq. (15). If we limit our deformations to terms which 
are quadratic in the Majorana operators, we can derive a general expression for the zero edge mode in the first-order 
approximation. To begin, let us express the dimer model (Hq) and the deformation (V) as: 

i i 

Hq = ^ ~ 2 '^i^ijlj^ (^ 3 ) 

with A and B real, skew-symmetric, and even-dimensional matrices. It follows that the eigenvalues of A (and B) are 
completely imaginary, and come in complex-conjugate pairs. The eigenproblem [i7o,7j] = thus has 

real eigenvalues Cn- [^0, '^n] = 2Cn'^n- The eigenvectors Vn form a complete basis for all linear Majorana operators, thus 
we may define the expansion 7^ = J]]/ ZjiVi] alternatively, Zji may be obtained from diagonalizing A. The eigenvalues 
{Cj} of —iA are interpreted as half the single-particle energies of Hq, i.e., Hq has single-particle excitations with energy 
2(j. Since Hq is the topological dimer model described in Eq. (15), we know that A has two zero modes (denoted 
Cl = C 2 =0)5 corresponding to the edge operators vi = 71 and V 2 = 72L- All other single-particle excitations are 
gapped, i.e., |2Cj| > Fq for j > 2, where Fq is the many-body spectral gap of Hq. The edge mode is derived as 

= - 2i y; y; BimZmnVn / = -2irQ y; y; Blm^mn1'n-F(2Cn/ro). (D4) 

m7^1,2L n7^1,2 m^l,2L n^l,2 


Since V is skew-symmetric, the sum over m excludes I; this sum also excludes 2L, since j 2 L is a zero mode and F 
an odd function (recall Sec. IV). Since Bq does not contain a term proportional to either 71 or 72L, it follows that 
the expansion jrn^i, 2 L = 2 excludes the zero modes vi = 71 and V 2 = 72L- Alternatively stated, all the 

relevant frequencies are gapped, i.e., |2Cn/ro| > I- Now applying that = —1/0. for these relevant frequencies, 

we find 

y; BimZmn/^- (D5) 

Applying this to the quadratic Majorana model (15), we obtain i[T>07 7i] = 73? as we have previously derived in 
less-general fashion; cf. Sec. IV. 


2. Deformed zero edge mode for the Z 3 chiral parafermion 


For illustration, we consider the Z3 chiral parafermion with 0 < 0 < tt/S. Recall that the filter function has not 
been specified in Eq. (C7). If we choose F{Cj) = we alternately derive Eendley’s edge mode. Substituting this 

choice into Eq. (C7) and applying Eq. (C6), we obtain 


— I V-> Vn I/I I I\ /I (jJ /I (jJ \ ii 


(D6) 


Here, we have introduced X (to distinguish Eendley’s choice) in place of X (which we reserve for the conventional 
quasi-adiabatic continuation). A useful relation is C 1 C 2 C 3 = 2 sin (30)/(3v^), which is derived from taking the 
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determinant of A4 as defined in Eq. (C5). This relation, in combination with elementary trigonometric identities, 
leads to 

(72 + (D7) 

The expressions for X and Y = —X\ in combination with Eq. (C 2 ) and (20), lead to the final result 


Oi =ji - 


IS 


2 sin (30) [ 


72 + A^73 + A + W 7 I ( 7 I + A ^ 7 ] + ^^^ 7273 ) 


(D 8 ) 


This deformed edge mode Oi is identical to the expression alternatively derived in Ref. 2, up to minor typographic 
errors in that reference. At least to first order in 5 , Fendley’s edge mode can thus be understood as a deformed quasi- 
adiabatic continuation of 71 ; the resultant operator is not optimally localized. This deformation cannot always be 
justified: as 0 (and Co) tends to zero, the localization of Oi worsens dramatically, as evidenced by the diverging 
first-order coefficient in Eq. (D 8 ), and the increasingly singular F. In contrast, the conventional quasi-adiabatic 
continuation chooses a filter function that is everywhere infinitely differentiable;^^ the resultant edge mode is well- 
defined for any 0 , as shown in App. C. 


Appendix E: Closing the chain by imposing translational invariance 

Our goal is to map a parafermion Hamiltonian (i^op) 011 an open chain of L sites, to a parafermion Hamiltonian 
(/C*^"^) on a closed ring. For simplicity, we assume that the open-chain Hamiltonian is translational-invariant up to 
edge corrections, i.e., we can decompose = X]r>i ^j,r^ where r is the range of the operator Vj^r^ i-e., V^-,r 

acts nontrivially on sites {j, j -h 1,..., j + r — 1}. Furthermore, we assume Vj^r is related to by translation. For 

a Z 3 example, consider H^p = — <^ 72 j+i 72 j + h.c., where Vj ^2 = —^ 72 j+i 72 i + h.c.. For N > 2^ a, Hamiltonian 

on a closed chain cannot be uniquely defined by the identification j 2 L+j = 7^5 due to the nontrivial commutation ( 1 ).^^ 
In our example, we might have tried to extend H^p by translational symmetry: H^p H^p + Vl ,2 5 then identified 
72 L +1 = 7 i directly. The result would be a closed-chain Hamiltonian 

L-l 

= - j y; + h.c. (El) 

i=i 

with a = 2 . On the other hand, first interchanging 7L+172E = ^72e72l+15 and then making the same identification, we 
arrive at a different Hamiltonian (El) with a = 1. If there exist more than one inter-edge couplings, the phase of each 
is ambiguous. These ambiguities motivate us to identify a^+j = cFj instead - since clock operators on different sites 
commute, the clock Hamiltonian is uniquely defined on a closed chain. Then by attachment of Jordan-Wigner strings 
we convert this clock Hamiltonian (which is clock-local) to a parafermion Hamiltonian (which is parafermion-local 
but clock-nonlocal). Let us illustrate this procedure with the family of open-chain Hamiltonians (from Eq. ( 2 )), 

which we extend and then identify = cti, to obtain clocks on a ring: 

N-l 

(E2) 

/3=0 

In the next step, we decompose the inter-edge coupling into a form such that Bi ( resp. B^) is a clock-local 

operator with a definite charge, and with support near site I (resp. L). To convert B^i to an inter-edge coupling (Hpf) 
of parafermions, we must attach as many powers of the string Q as the charge of B^. The resultant inter-edge coupling 

^7 =B;B,(w“g)[[®^ll (E3) 

can always be expressed locally in terms of parafermions, as follows from the Jordan-Wigner transformation (5). Here 
in Eq. (E3), [[Br]] refers to the charge of B^^ i.e., B^Q = and we have also introduced the twist parameter 

a G Z]sf. In our example (E2), [[Br]] = [[(Tl]] = I, so we attach aL ctlQ ^ 72 l- Then from ctlctiQ = 721.5 

we obtain parafermions-on-a-ring: 

N-l 

/C7 = - J ^ ■ (E4) 

/3=0 
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Appendix F: Deriving the topological order parameter in the XYZ model 

Our aim is to the derive the topological order parameter for the XYZ Majorana model (28), which we parametrize 
by Jx and Jy. There are two parameter regions with distinct order parameters: (i) | Ja,|, | < 1, as described in App. 

FI, and (ii) \Jy \ > 1, \ Jy \ > \ Jx\, described in App. F2. 


1. Topological order parameter for | Jcc|, | Jy| < 1 

Let us parametrize the XYZ model (28) as 


v-A J 

Hs = Ho + sV, where F = X ('T23-i72,-72i+i72j+2 - «i72^-i72i+2 ), s = Jxt=^, (FI) 

i=i “ 

and the frustration-free Hq is defined in Eq. (26). Hq commutes with the zero edge modes 7i and 721,, thus its 
groundstate space is described by a topological order parameter Q = —271721.. If we frustrate this model through E, 
the spectral gap above the lowest two states remain finite in the parameter region: | < 1- Our task is to derive 

the dressed order parameter in this region. We will employ the identities: 

[r7i] =-2727374 + 2it74 = (-7273 + + u'-), and [V, 72^] = (2++ + 72^-272^-1), (F2) 

where w± = 74 ± and z± = 72i_3 ± i'yiL-i are eigenoperators under commutation by Hq, i.e., [Ho,w±] = ^2w±, 
and [Hq,z±] = ±2z±. It follows that 

*[F^o, 7 i] = y = {it - 7273)(i75) and 

*[^*0,721,] = J dt-F(2t)e®-^°*[y,72i,]e“*-^“* = -n2i,-4(-*^ + 721,-2721,-1), (F3) 

which determine the edge modes Oi and Or to first order in s; see Eq. (20). Applying that Q = —iOiOr, we are 
immediately led to Eq. (30). In this derivation, we have assumed that \Jx\ < 1, Tc 7^ 0; if instead Jx = 0 and 
\^y I < 1, Jy 0^ the final answer (30) is identical. 


2. Deriving the topological order parameter for \Jy I ^ Ijd^yl > \Jx\ 

Up to a proportionality constant, we parametrize the XYZ model (28) as iL' = Hq + sU', where s = I/Jy, 


L-l L-l L-l 

Hq = and V = J:r^l2j-il2jl^j+il2o+2 + i^l202j+i- (F4) 

j = l j = l j = l 

The frustration-free commutes with the zero edge modes 72 and 721.-1, thus its groundstate space is described by a 
topological order parameter = —i^ 2 l 2 L-i^ which must be distinguished from the order parameter (Q = of 

iLo- If we frustrate this model through U', the spectral gap above the lowest two states remain finite in the parameter 
region: \Jy I > ^My\ > \Jx\- In this region, the order parameter is dressed as 

Jx ^ 

Qy = -^727212-1 + ^(7174767211-1 + 727211-57211-37211) + -r{lel2L-i + 727211-5) + • • • (F5) 

Jy Jy 

Its derivation straightforwardly generalizes App. El. Defining as the first-order truncation of (i.e., dropping 
the dots), we define a closed-chain Hamiltonian by Hxyz — { — where Hxyz is defined in Eq. (28). We 

numerically evaluate the groundstate parities of this closed-chain Hamiltonian for both a G Z2, on a chain of 22 sites. 
The result is shown in Eig. 4(d); we have whitened the parameter regions where the groundstate parity switches upon 
fiux insertion. Note that \Jy \ > I, | > \Jx\ corresponds to regions B and D in this figure. 
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Appendix G: Finite-size analysis of the translational-invariant Majorana chain 


In this Appendix, we present a finite-size analysis of the XYZ Majorana model (28) for the parameter region: 
\Jx\') \Jy I < 1. which we label as the square C in Fig. 4. Suppose we close the chain with the translational-invariant 
boundary conditions (TIBC); a detailed description of this procedure is available in App. E. We then obtain two 
Hamiltonians, which we denote by ICxyz (resp. ICxyz) in the case of periodic (resp. antiperiodic) boundary condition. 
A signature of topological order is that the groundstate of ICxyz is sensitive to its boundary condition; in particular, 
we expect that the groundstate parity switches between ICxyz and ICxyz^ for all parameters within the square C. A 
finite-size study shows that the parity is invariant for a subregion (colored black) of C where quantum fluctuations 
are dominant; see Fig. 4(b). The aim of this Appendix is to explain the effect of finite size on the groundstate 
sensitivity. Our explanation relies on a Majorana-spin duality which we will first introduce; one motivation for the 
spin representation is that numerical simulations are often easier with a clock-local basis. 


(i) To describe this duality, we first relate the open-chain Majorana Hamiltonian {Hxyz) to a closed-chain spin 
Hamiltonian (Qxyz) fo App. G 1. We argue that Hxyz being topologically ordered implies clock-local order in Qxyz- 

(ii) Then in App. G2, we show that the closed-chain Majorana Hamiltonian (IC^xyz) is dual to the closed-chain spin 
Hamiltonian (Gxyz)- Applying this duality, we explain what clock-local order in Qxyz implies for IC^xlz- We then 
analyze the groundstates of IC^xlz and their sensitivity to twisting. 


(iii) Finally in App. G 3, we describe a finite-size comparative study of the three types of boundary conditions discussed 
in this paper: (i) TIBC, (ii) zeroth-order topological boundary condition (TBC), and (iii) first-order TBC. 


1 . Relating topological order in Hxyz to clock-local order in the spin Hamiltonian Qxyz 
As defined in Eq. (28), the open-chain Majorana Hamiltonian can be rewritten as the Heisenberg spin model: 

L-l 

H^rz = -Yl E with J, = 1. (Gl) 


Here, are Pauli matrices acting on site j. To derive the above spin representation, we identify = a, = r, and 
apply the Jordan-Wigner transformation (5): for bulk couplings (1 < j < L), 


44+1 =-T'2.-i72.%.+i 72,+2> and = (G2) 

We define Qxyz as the translational-invariant extension of Hxyz to a ring: 

Qxyz = ^ (G3) 


Here, periodic spin boundary conditions are imposed as: = 4^ The phase diagram of Qxyz is well-understood 

analytically. Qxyz is gapless where | J^| = \Jjy\ > | J^|, i.e., for any two parameters having equal magnitude that 
is greater than or equal to the third magnitude; otherwise, the phase is gapped with a clock-local order parameter 4^, 
corresponding to of the largest magnitude. We illustrate this symmetry-breaking in Fig. 8(a), for the parameters 

(J., = (0-1.-0-2,1). 
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(a)//^, (h)g^, (c)ICZ (d)/CSz 



(^^^HxyZ 


mg^z (g)/Cz wjc^z 



FIG. 8. Lowest-lying spectra of the XYZ model on a chain of 10 sites. Two parametrizations are considered: (a)-(d) cor¬ 
respond to the same parameters [Jx^ Jy^ Jz) — (0.1,—0.2,1), while (e)-(h) correspond to (Jx, Jy, Jz) = (0.75,-0.75,1). The 
spectrum of each Hamiltonian is split into even (blue) and odd (red) charge sectors, (a) and (e): spectrum of the open-chain 
HxYz{Jx,Jy,Jz), (b) and (f): Qxyz, ( c ) and (g): IC^xyz^ (d) and (h): IC^xyz- 

We argue that Hxyz being topologically ordered (in the fermionic language) indicates that Gxyz is symmetry- 
broken (in the spin language). This implies that Hxyz and Gxyz simultaneously manifest a degenerate groundstate 
manifold, as illustrated in Fig. 8(a) and (b) for the parameters ( Jy^Jz) = (0.1, —0.2,1). This follows because Gxyz 
is the translational-invariant extension of Hxyzi thus both Hamiltonians are approximately minimized by the same 
wavefunctions. Alternatively stated, we expect that the spectral gaps of Hxyz and Gxyz (above their two lowest-lying 
states) are related multiplicatively - a quantum phase transition in Hxyz indicates a simultaneous transition in Gxyz- 
To support our hypothesis, we numerically evaluate the spectral gap of Hxyz in Fig- 4(a); the gapless lines (colored 
black) separate the plot into five regions labelled by A to with 5, C and D corresponding to topological phases 
(cf. Sec. V). Our hypothesis is supported by the exact coincidence of gapless lines in both Hamiltonians {Hxyz and 
Gxyz)- Here, we are comparing our numerical simulation of Hxyz with a known analytical result about Gxyz] cf. the 
previous paragraph. 


2. Majorana-spin duality on a closed chain 

We extend the open-chain Hxyz onto a ring through the prescription (E3), to obtain the closed-chain Hamiltonian 
IG^xyz = Hxyz + ^H^xyz^ as shown in Eq. (35). The inter-edge coupling can be rewritten as 

- Jy^l4{-irQ - Jxi^li-ITQ. (G4) 

Let us define for any Hamiltonian H a projected Hamiltonian [H]f 3 in the charge sector i.e., [H]f 3 acts in the space 
of functions with eigenvalue under conjugation by Q. Clearly, 

(-1)“+^. (G5) 

By inspecting the matrix elements of Hxyz^ Gxyz and IC^xlz ia each charge sector, we derive a duality between 
Majorana fermions and spins on a closed chain: 

[fG^^XYzja = [GxYz]a^ for a G Z2. (G6) 

This duality on a closed chain generalizes an open-chain duality^^’^^ between the transverse-field Ising model and the 

Kitaev wire. Gxyz being translational-invariant then also implies the same for ICxyz (resp. ICxyz) in the even (resp. 
odd) charge sector, as we alluded to in Sec. VB. 

This duality is further illustrated in Eig. 8(b)-(d) for the parameters {Jxi Jy^ Jz) = (0.1,—0.2,1). If we define 
as the groundstate of Gxyz in the charge sector a G Z2, then clearly 

minimizes [K}xyz]o, and minimizes [K}xyz]i- (G7) 
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For ( Jy, J^) = (0.1, -0.2,1), 


minimizes fC^xlz^ and also minimizes 1C 


( 1 ) 


(G8) 


This is illustrated in the spectra of Fig. 8(b)-(d), and explains why the groundstate parity switches as we twist 1C 


(a) 

XYZ - 


Due to the duality (G6), the first statement (G7) is exactly satisfied for any finite size; the second statement 
(G8) is not necessarily satisfied on finite-size chains, i.e., the actual groundstate of IC^xlz need not originate from 
the symmetry-broken manifold of Qxyz- For example, we consider the spectra at ( Jy, Jz) = (0.75, —0.75,1) for a 
10-site chain, as plotted in Fig. 8(e),(f),(g) and (h). Fig. 8(f) demonstrates the finite-size splitting of the low-energy 
subspace of Qxyz] in particular, i^ energetically favored in comparison with While 

also minimizes ICxyz (see Fig. 8(g)), does not minimize Kxyz (see Fig. 8(h)). A corollary is that the 

parity does not switch, and we have verified that this behavior persists for a chain of 30 sites. 


As a side remark, we note that the duality described in Eq. (G6) is easily generalized to a parafermion-clock duality: 
if IC^"^^ is the Hamiltonian for parafermions on a ring, with a particular twist a G Z^v, and Q is the Hamiltonian for 
clocks on a ring, then 


[IC^^ = [Q]a, with a G Zat. 


(G9) 


Let us exemplify this duality for the frustration-free models i7(iv,i) defined in Eq. (2). The closed-chain parafermion 
Hamiltonian (E4) is dual to the closed-chain clock Hamiltonian (0(iv,i)) Eq. (E2). That is, 1C^n\ and Q^n,i) satisfy 
Eq. (G9), as follows from the identity (25). 


3. Comparing the three types of closed-chain boundary conditions 


Applying the parity switch as a criterion, we would like to evaluate the various types of closed-chain boundary 
conditions. We focus on the Jx = —Jy line, for Jx G [0, 1), and define Jx such that the parity switches for 0 < 
and is invariant for Jx < Jx < With the topological boundary condition (TBG) of Eq. (32), the parity switches 
for all parameters where the open-chain Hamiltonian is topological, i.e., Jx = 1, as proven in Sec. V A. With the 
translational-invariant boundary condition (TIBG) of Eq. (E3), it has also been argued^^ that Jx = ^ for system sizes 
much larger than the correlation length; in practice this can only be verified with a much larger system size than can 
be simulated with Lanzcos. In Eig. 9, we numerically evaluate Jx for three different types of boundary conditions: (i) 
TIBG, (ii) zeroth-order TBG (see Sec. V A), and (hi) first-order TBG (see Eq. (33)). Eor all the system sizes that we 
simulate (up to L = 30), the groundstates with (i) and (ii) are comparably sensitive to twisting, while the groundstate 
with (iii) is most sensitive. A naive extrapolation suggests that these conclusions are robust in the thermodynamic 
limit. 
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FIG. 9. Finite-size analysis of the closed-chain XYZ model along the Jx — —Jy line, for Jx G [0,1). We close the chain using 
three types of boundary conditions: (i) translational-invariant (colored red), (ii) first-order topological (blue), and (hi) zeroth- 
order topological (green). We define Jx such that upon twisting the boundary conditions, the groundstate parity switches for 
0 < Jcc < but is invariant for < 1. We numerically evaluate Jx for varying system sizes (up to L = 30), as shown by 

solid dots in the plot. Each dashed line is a fit to the function: a + 6/L + c/L^, and by extrapolation we obtain an estimate of Jx 
in the thermodynamic limit. Since the fitting function is not theoretically motivated, we emphasize that these extrapolations 
are more suggestive than rigorous. 


Appendix H: Model with coexisting topological and symmetry-breaking orders 

We present a model example with coexisting topological and symmetry-breaking orders; all other examples in this 
paper have either been purely-topological or purely-symmetry-breaking. We consider Zg parafermions which have 
four distinct phases labelled by the possible integer divisors: n G {1,2,4, 8}; each of these phases may be realized 
by the frustration-free models in Eq. (37). Defining m = S/n and g =gcd(m,n), we consider the phase labelled by 
(n = 2, m = 4, ^ = 2). By the general classification proposed in Ref. 30, > 1 implies there is symmetry-breaking, 

and m > g implies topological order; the rest of this Appendix is an elaboration of these statements, with emphasis 
on concepts we have introduced in this paper, namely their topological edge modes and their local indistinguishability. 

We begin with the most intuitive groundstate basis of 77(8,2) in Eq. (37) - these are fully polarized as 

|'0a,2) = \oi). + |(a + 4 ). ), (HI) 

i=i ^ 

with Of G Z4 = {0,1,2,3} labelling four degenerate groundstates on an open chain. Recall in the above expression 
that cFj\a)^ = A simple computation shows that and are clock-local order parameters for any subcell 

coordinate j : 


Tf\^Pa,2) = \^a,2) and <j]\tpa,2) = l^'^°‘\'tpa,2) ■ (H2) 

This implies the parafermionic order parameters for any j, since J 2 j-i ^ ^j^j-i^j -2 • • • J 2 j ^ ^j^j^j-i^j -2 • • - ^ 
through the Jordan-Wigner transformation. While in this model ^ = 2 > 1, we might ask more generally why 
^ > 1 is necessary for the existence of a parafermionic order parameter. This answer is found in Ref. 30, as we 
summarize: for the family of Hamiltonians parametrized by {N^n) in Eq. (37), the generalization of Eq. (H2) is that 
and are clock-local order parameters - then if there exist integers a and b such that am = bn but is not a 
multiple of A, then one might string together the following order parameter - ^ which up to a phase 

factor may be expressed as a parafermionic order parameter: with c G {0,1,..., ^ —1}; clearly here we need g > 1. 

Returning to our model example, we would like to show that this four-fold degeneracy is really a product of two 
degeneracies which each arises from a different kind of order - the topological degeneracy is attributed to edge modes, 
while the symmetry-broken degeneracy is attributed to the parafermionic order parameter. One way to proceed 
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is to apply the parafermionic order parameter as a symmetry-breaking ‘field’ - then if there exists any remnant 
degeneracy, it must arise from topological order. As an intermediate step in this computation, we introduce a basis 
that diagonalizes the Zg generator (recall the expression for Q in Eq. (3)): 

= (H3) 

/3=0 

which satisfies 

Q|</'a,2) = W^“|(/)a,2), r/|(/)a,2) = |0a,2) and CT||(/)a,2) = |0a-lmod4,2)- (H4) 

The effect of the symmetry-breaking ‘field’ is to resolve the groundstate manifold as: 

\O 0 ,s) = T( + {-lf\cl>s+ 2 , 2 } ) (H5) 


with /3, ^ G Z 2 ; this basis splits as 


= (H 6 ) 

In this sense, is a symmetry-breaking index that is distinguished by the parafermionic order parameter, while S is 
a topological index that labels each state in each symmetry-broken manifold. While here ^ G Z 2 , in the most general 
case 6 G Z^/^ and therefore there is no remnant topological degeneracy unless m > g. Returning to our model, while 
7 ^ breaks the Zg symmetry generated by Q, there remains a Z4 symmetry generated by and a simple computation 
shows that \0p^s) is an eigenstate of with eigenvalue (—1)'^. This remnant symmetry generator has the following 
fractionalized representation in the groundstate space: 

Q 2 = ^7r7lL = (H7) 

This operator has the same groundstate matrix elements as as follows from acting trivially in the ground- 

state space. Indeed, from Eq. (H4), crG^cr|| 0 cK, 2 ) = | 0 a, 2 ) for any i, j. Within each symmetry-broken manifold labelled 
by yd, the topological edge modes 7 ^ and 7 !^ permute the groundstates as (/d, ^ + lmod 2 ). Alternatively 

stated, each symmetry-broken manifold forms a representation of the non-commutative algebra: 

7?72L = - 7|l7i , 7i Q 2 = - Q 271 , 

and 7 |iQ 2 = -Q 272 L- (H 8 ) 

While the discussion thus far applies to an open chain, we observe that Q 2 is a legal Hamiltonian term on a closed 
chain, and its application will split this topological degeneracy. 

It is instructive to consider a closed-chain Hamiltonian without both symmetry-breaking ‘field’ and topological edge 
modes. Since the ‘field’ is an operator with quartic charge (i.e., and | 0 a, 2 ) has charge 2a (cf. Eq. 

(H4)), we have that the ‘field’ permutes the groundstates as \4>a,2) | 0 a- 2 , 2 )- The open-chain manifold then divides 

naturally into two subspaces which are individually condensates of the ‘field’: 

Ps = \4’S,2){4>6,2\ + \4>2+6,2){4>2+S,2\, with S G Z 2 . (H9) 

This reveals another perspective to the two different kinds of degeneracies: the topological (symmetry-broken) de¬ 
generacy depends (is independent of) the spatial topology; the net degeneracy on an open chain is a product of both 
degeneracies. 


Appendix I: Indistinguishability by symmetric, parafermion-local operators 

Let |0aez^,n) denote the groundstates of the frustration-free Hamiltonian i7(iv,n)- Our basis is chosen as Q|0a,n) = 
Assuming g = gcd(m,n) > 1, divide P[i 7 ( 7 v,n)] into (m/g) subspaces labelled by ^ G Z^/^: Ps = 
\^am/g-\-S,n){4^am/g-\-S,n\‘ Let US usc O to denote a Symmetric, parafermion-local probe. We would like to show 
that the g x g matrix = (0 am/g-\-6,n m/g-\-6,n) is proportioual to the identity, up to superpolynomially-small 
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finite-size corrections; this is equivalent to the condition (49). The off-diagonal elements vanish since [0,Q] = 0. The 
difference between diagonal elements is 

^ m—1 

where |'0^ez^,ri) are fully-polarized states (39) related to |0) by the Fourier transform (40). Crucially, the differ¬ 
ence in phase factors vanish when fi = v mod thus M.^aa — reduces to a sum of terms proportional to 

7^,zy = (' 0 /x,n| 0 |' 0 zy,n )5 with ji ^ h' mod Q. A Straightforward generalization of App. B shows that vanishes 
superpolynomially in the system size, if O is clock-local. Now we address the case when O is clock-nonlocal but 
parafermion-local, as is possible if O is an inter-edge coupling with support near sites 1 and L. Any inter-edge cou¬ 
pling may be decomposed into the form: where Bi (resp. Br) is a clock-local operator with a definite 

charge, and with support near site 1 (resp. L). The nonlocal clock-parafermion transformation (5) implies that there 
are as many powers of the string Q as the charge of B^. Once again, our notation for the charge is: B^Q = 

We provide an example of this decomposition for Z 4 parafermions: 


I1I2I2L-1J2L = wcTiTicrlriQ^, ( 12 ) 

where Bl = ujcfiTi^ Br = cf\tl and [[Br]] = 2 . 

We then claim that 


\'lpu,n) = 0, for £ Zm, I' Hiod g, (13) 

and for finite-ranged Br and Bi; our claim will be proven in App. II. Eq. (13) implies the vanishing of Eq. (II), 
and thus a proof of Eq. (49) in the frustration-free limit. If Bi and Br are superpolynomially decaying into the bulk, 
Eq. (II) vanishes to up to superpolynomially-small finite-size corrections, which can be bound in a manner analogous 
to App. B. Einally, let us consider deviating symmetrically from the frustration-free limit, while preserving the gap 
within each of the g charge sectors of Ps. Then the deformed subspace is related to the old by a quasi-adiabatic 
continuation which is distinct for each ^-multiplet labelled by 6: Ps^s = Then applying [Vs^siQ] = 0? 


Ps,sOPs,s = Vs,sPsVldVs,sPsVl^ = Vs,sPs{VlBiVs,s){VlBrVs,s)Q^^^^^'^PsV^ 


t 

6,s' 


(14) 


Since Vs^s is locality-preserving, the dressed operators ^BiVs^s and ^BrVs^s are superpolynomially-decaying into 

the bulk, thus we may apply our frustration-free result: Ps{V^^BiVs^s){ys ^ final result is 

that the indistinguishability condition is robust to superpolynomial accuracy. 


1. Vanishing of Eq. (13) for finite-ranged Bi and Br 


Let us decompose BiBrQ^^^'^^'^ into a sum of terms like Oj, where r is a nonzero complex number, and Oj 

is an on-site operator acting on site j, with the general form for integers aj and bj. In our Z 4 example (12), 


cr|r| ^ aj = 2, bj = 3 

Tj ^ aj = 0, bj = 2 


for jG{l,L} 

for j G {2,3,...,L-2,L- 1} 


(15) 


We deduce two useful properties of {a^-, bj} from the general form of BiBrQ^^^'^'^^: 

(a) Br has support over a set of sites near L, thus we may express Br oc n^=i.-fi+i some range (R) that is 

much smaller than the length of the chain. Since [[(Jj]] = 1 while [[rj]] = 0 , [[Br]] = example, 

Br = and [[Br]] = aL = 2. 


(b) Bl and Br has no support on any site j far away from any edge, where the only nontrivial operation 
arises due to the string operator Thus for j in the bulk, bj = [[Br]] and aj = 0. In example (15), the bulk 

sites are labelled hy j G {2, 3,..., L — 2, L — 1}. 
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Applying Eq. (39), we now express Eq. (13) as a sum of products of on-site elements, with the form: 


L ^ n—1 

= where \a,n). = + mP)aj\a) . = oj°‘\a) 

3 = 1 /3=0 


(16) 


Each on-site element may be evaluated as 


lu. 


’’A= E 

Ck:,/3= 
n—1 

= E 




-(/i + ma V -\-bj m^') 


Ck:,/3=0 

q :=0 


^ (/x+ma) 


mod m — ^ ^/x,zx+6j mod m ^aj,0modr 


(17) 


The second equality follows from /i and z/ G In our Z4 example, where g = m = n = 2, 

(0,2|Oi|i,2). 

^ f5((0|i + + |3)i) = 5((0|i + (2|j)(|0)j + |2)j) = 1, for je{l,L} 

U((0|, + (2|pry|l)^, + |3)p = i((0|^, + (2|p(|l). + |3)p = 0, for j e {2,3,..., L - 1}. 

Erom the vanishing of on-site elements in the bulk, we deduce that 7oi = 0. Suppose we try to have the bulk terms 
be non-vanishing; in the following we show that one of the edge terms must consequently vanish! Erom the second 
line of Eq. (18), we would need a bulk operator with odd bj. Erom property (b) above, we would then need that 
[[Br]] is also odd. This implies from property (a) above that ^ range of the Hamiltonian, 

i.e., for at least one site labelled by j G {L — + 1,..., L}, aj is odd. Eor this particular site, we evaluate 

(0,2|qHt2). = i((0|. + (2|.)afr^d|l)j + |3)j). (19) 

If bj is even, then this term reduces to a sum of (0|l), (0|3), (2|l), and (2|3), which individually vanishes. Suppose 
we choose bj to be odd, then 

{0,2\Oj\i,2)-. = 1((0|. + (2|p<Tf (|0). + \2)-.) = i((0|. + (2|.)(|0)j + (-1)“^12).) = 0. (110) 

The last equality follows from aj being odd. We have thus proven that one edge term must vanish, and therefore 
7oi = 0 most generally. 


More generally, we would like to show that = 0 for /r 7^ z/ mod g. If we assume the converse, we require that 
each on-site element must be non-vanishing, i.e., from Eq. (17), 

\fj G {1,2,..., I/}, g = u ^bj mod m, and aj = 0 mod n. (HI) 

Erom our assumption, we now prove a contradiction: 

(i) Eor a site j far away from the edge, we derive from Eq. (Ill) and property (b) that [[Br]] = g — v mod m. The 
condition on g and z/ then implies that [[Br]] is not an integer multiple of g. 

(ii) / 0 implies from Eq. (17) that each of Oj is a multiple of n. With property (a), we derive that [[Br]] = 
'^^=L-R+i % ~ where k G Z^. Since g = gcd(m,n), this implies [[Br]] is an integer multiple of g, in contradiction 
with (i). Thus we conclude that = 0 for /r 7^ z/ mod g. 


^ G. Moore and N. Read, Nucl. Phys. B 360 , 362 (1991). 
^ N. Read and D. Green, Phys. Rev. B 61 , 10267 (2000). 


^ E. Majorana, Nuovo Gimento 5, 171 (1937). 
2 P. Fendley, J. Stat. Mech. , P11020 (2012). 



27 


^ D. A. Ivanov, Phys. Rev. Lett. 86, 268 (2001). 

® A. Y. Kitaev, Phys.-Usp. 44, 131 (2001). 

^ A. Y. Kitaev, Ann. Phys. 1, 2 (2003). 

® C. Nayak, S. H. Simon, A. Stern, M. Freedman, and S. D. 

Sarma, Rev. Mod. Phys. 80, 1083 (2008). 

^ J. Alicea, Y. Oreg, G. Refael, F. von Oppen, and M. P. A. 
Fisher, Nature Phys. 7, 412 (2011). 

S. Ostlund, Phys. Rev. B 24, 398 (1981). 

“ D. A. Huse, Phys. Rev. B 24, 5180 (1981). 

D. A. Huse, A. M. Szpilka, and M. E. Fisher, Physica A 
121, 363 (1983). 

L. K. S. Howes and M. den Nijs, Nucl. Phys. B 215, 169 
(1983). 

G. Albert ini, B. M. McGoy, and J. Perk, Phys. Lett. A 
139, 204 (1989). 

H. An-Yang and J. Perk, Intern. J. Mod. Phys. B 11, 11 
(1997). 

R. B. Potts, Proc. Gamb. Phil. Soc. 48, 10 (1952). 

S. Elitzur, R. B. Pearson, and J. Shigemitsu, Phys. Rev. 
D 19, 3698 (1979). 

G. W. Semenoff and P. Sodano, J. Phys. B: At. Mol. Opt. 
Phys. 40, 1479 (2007). 

L. Eu and G. L. Kane, Phys. Rev. Lett. 100, 096407 (2008); 
Phys. Rev. B 79, 161408 (2009). 

R. M. Lutchyn, J. D. San, and S. D. Sarma, Phys. Rev. 
Lett. 105, 077001 (2010). 

Y. Oreg, G. Refael, and E. von Oppen, Phys. Rev. Lett. 
105, 177002 (2010). 

A. G. Potter and P. A. Lee, Phys. Rev. Lett. 105, 227003 

( 2010 ). 

T. D. Stanescu and S. Tewari, J. Phys.: Gondens. Matter 
25, 233201 (2013). 

V. Mourik, K. Zuo, S. Erolov, S. Plissard, E. Bakkers, and 
L. Kouwenhoven, Science 336, 1003 (2012). 

A. Das, Y. Rouen, Y. Most, Y. Oreg, M. Heiblum, and 

H. Shtrikman, Nature Phys. 8, 887 (2012). 

S. Nadj-Perge, 1. K. Drozdov, J. Li, H. Ghen, S. Jeon, 
J. Seo, A. H. MacDonald, B. A. Bernevig, and A. Yazdani, 
Science 346, 602 (2014). 

A. M. Turner, E. Pollmann, and E. Berg, Phys. Rev. B 
83, 075102 (2011). 

J. Motruk, A. M. Turner, E. Berg, and E. Pollmann, Phys. 
Rev. B 88, 085115 (2013). 

G. Kells, arXiv: 1409.6575. 

R. Bondesan and T. Quella, J. Stat. Mech. , P10024 (2013). 

L. Eidkowski, R. M. Lutchyn, G. Nayak, and M. P. Eisher, 
Phys. Rev. B 84, 195436 (2011). 

M. P. Zaletel, R. S. K. Mong, and E. Pollmann, Journal 
of Statistical Mechanics: Theory and Experiment 2014, 
P10007 (2014). 


R. B. Laughlin, Phys. Rev. B 23, 5632 (1981). 

X. Ghen, Z.-G. Gu, and X.-G. Wen, Phys. Rev. B 83, 
035107 (2011). 

E. Eradkin and L. P. Kadanoff, Nucl. Phys. B 170, 1 
(1980). 

B. Nachtergaele, Gommun. Math. Phys. 175, 565 (1996). 
A. S. Jermyn, R. S. K. Mong, J. Alicea, and P. Eendley, 
Phys. Rev. B 90, 165106 (2014). 

M. B. Hastings and X.-G. Wen, Phys. Rev. B 72, 045141 
(2005). 

M. B. Hastings, arXiv:1001.5280. 

S. Bravyi and M. Hastings, Gommun. Math. Phys. 307, 
609 (2011). 

H.-J. Mikeska and A. K. Kolezhuk, Lecture Notes in 
Physics 645, 1 (2004). 

A. Kitaev and J. Preskill, Phys. Rev. Lett. 96, 110404 
(2006). 

M. Levin and X.-G. Wen, Phys. Rev. Lett. 96, 110405 
(2006). 

H. Li and E. D. M. Haldane, Phys. Rev. Lett. 101, 010504 
(2008). 

A. Alexandradinata, T. L. Hughes, and B. A. Bernevig, 
Phys. Rev. B 84, 195103 (2011). 

D. J. Glarke, J. Alicea, and K. Shtengel, Nature Gommu- 
nications 4, 1348 (2011). 

N. H. Lindner, E. Berg, G. Refael, and A. Stern, Phys. 
Rev. X 2, 041002 (2012). 

M. Gheng, Phys. Rev. B 86, 195126 (2012). 

M. B. Hastings, arXiv:1008.5137. 

S. Bravyi, M. Hastings, and S. Michalakis, J. Math. Phys. 
51, 093512 (2010). 

M. A. Levin and X.-G. Wen, Phys. Rev. B 71, 045110 
(2005). 

S. Michalakis and J. Pytel, Gommun. Math. Phys. 322, 
277 (2013). 

G. Ortiz, J. Dukelsky, E. Gobanera, G. Esebbag, and 
G. Beenakker, Phys. Rev. Lett. 113, 267002 (2014). 

Y. Zhuang, H. J. Ghanglani, N. M. Tubman, and 

T. Hughes, arXiv: 1502.05049. 

E. H. Lieb, Gommun. Math. Phys. 28, 251 (1972). 

A. E. Ingham, J. London Math. Soc. 9, 29 (1934). 

L. A. Takhtadzhan and L. D. Eaddeev, Russ. Math. Surv. 
34, 11 (1979). 

E. Ercolessi, S. Evangelisti, E. Eranchini, and E. Ravanini, 
Phys. Rev. B 88, 104418 (2013). 

A. Kitaev and G. Laumann, arXiv:0904.2771. 

M. Greiter, V. Schnells, and R. Thomale, Annals of 
Physics 351, 1026 (2014). 



